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Abstract: We outline a general strategy for measuring spins, couplings and mixing angles 
in the case of a heavy partner decay chain terminating in an invisible particle. We consider 
the common example of a heavy scalar or fermion D decaying sequentially to other heavy 
particles C, B and A by emitting a quark jet j and two leptons and . We derive analytic 
formulas for the dilepton ({^^^~}) and the two jet-lepton ({j^n} and {j^/}) invariant mass 
distributions for the case of most general couplings and mixing angles of the heavy partners. 
We then consider various spin assignments for the heavy particles A, B, C and D, and 
for each case, derive the relevant functional basis for the invariant mass distributions which 
contains the intrinsic spin information and does not depend on the couplings and mixing 
angles. We propose a new method for determining the spins of the heavy partners, using the 
three experimentally observable distributions + and — {j^ }• We 

show that the former two only depend on a single model-dependent parameter a, while the 
latter may depend on two other parameters (3 and 7. By fitting these distributions to our set 
of basis functions, we are able to do a pure measurement of the spins per se. Our method 
is also applicable at a pp collider such as the Tevatron, for which the previously proposed 
lepton charge asymmetry is identically zero and does not contain any spin information. In 
the process of determining the spins, we also end up with an independent measurement of 
the parameters a, (3 and 7, which represent certain combinations of the couplings and the 
mixing angles of the heavy partners A, B, C and D. 
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1. Introduction 



The ongoing Run II of the Fermilab Tevatron and the imminent turn-on of the Large Hadron 
CoUider (LHC) at CERN are beginning to explore the physics of the Terascale. There are 
sound theoretical reasons to believe that some new physics beyond the Standard Model (BSM) 
is going to be revealed in those experiments. Perhaps the most compelling phenomenological 
evidence for BSM particles and interactions at the TeV scale is provided by the dark matter 
problem [1] . It is a tantalizing coincidence that a neutral, weakly interacting massive particle 
(WIMP) in the TeV range can explain all of the observed dark matter in the Universe. A 
typical WIMP does not interact in the detector and can only manifest itself as missing energy. 
The WIMP idea therefore greatly motivates the study of missing energy signatures at the 
Tevatron and the LHC [2]. 

The long lifetime of the dark matter WIMPs is typically ensured by some new exact 
symmetry, e.g. i?-parity in supersymmetry [3], KK parity in models with extra dimensions [4], 
T-parity in Little Higgs models [5,6], [/-parity [7,8] etc. The particles of the Standard Model 
(SM) are not charged under this new symmetry, but the new particles are, and the lightest 
among them is the dark matter WIMP. This setup guarantees that the WIMP cannot decay, 
and more importantly, that WIMPs are always pair-produced at colliders. The cross-sections 
for direct production of WIMPs (tagged with a jet or a photon from initial state radiation) at 
hadron colliders are typically too small to allow observation above the SM backgrounds [9]. 
Therefore one typically concentrates on the pair production of the other, heavier particles 
(e.g. superpartners, KK-partners, or T-partners), which also carry nontrivial new quantum 
numbers just like the WIMPs. Once produced, those heavier partners will cascade decay 
down, emitting SM particles which are in principle observable in the detector. However, each 
such cascade also inevitably ends up with an invisible WIMP, whose energy and momentum 
are unknown. Since the heavy partners are being pair-produced, there are two such cascades 
in each event, and therefore, two unknown WIMP momenta. In addition, at hadron colliders 
the total parton level energy and momentum in the center of mass frame are also unknown, 
and thus the exact reconstruction of the decay chains on an event by event basis is a very 
challenging task [10-12]. 

The lack of fully reconstructed events makes the mass and spin determination of the 
heavy partners rather difficult. Due to the escaping WIMPs, the heavy partners cannot be 
reconstructed as resonances in the invariant mass distributions of their decay products. Their 
masses therefore must be measured from (a sufficient number of) kinematic endpoints [13-17]. 
The method can be successful, if a suitable cascade decay chain is identified in the data. An 
example of such a decay chain is presented in Fig. [l|, where we show the sequence of three 
two-body decays D^C + q, C^B + £n and B — > A + if. Here D, C, B and A are 
some heavy particles with masses mo, mc, ms and uia, correspondingly. For simplicity, 
throughout this paper we shall assume that all heavy particles are on-shell, i.e. 

niD > rnc > ms > ruA ■ (1-1) 
We shall take the visible decay products to be a quark jet q and two leptons (either electron 
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Figure 1: The typical cascade decay chain under consideration in this paper. At each vertex we 
assume the most general coupling (see Sec. 1.2 for the exact definitions) and we quote our results in 
terms of the dimensionless mass ratios x, y and z. 



or muon), in that order^. For discussion purposes, the leptons are often referred to as "near" 
{in) and "far" (^/), although this distinction is difficult to make in the actual data. Our 
setup follows closely the conventions of Refs. [17,19-22]. Accordingly, we shall also find it 
convenient to express our results in terms of the mass ratios 

2 2 2 

x = -^, y = ^, z = ^. (1.2) 
mjj itIq 

For a variety of reasons, the particular decay sequence exhibited in Fig. |l| has attracted 
a lot of interest in the past and has been extensively studied both in relation to an even- 
tual discovery of new physics as well as precision measurements of the new physics param- 
eters. Rather early on, it was realized that this decay chain commonly occurs in the most 
popular models of low energy super symmetry, such as minimal supergravity (MSUGRA), 
minimal gauge mediation [23], minimal anomaly mediation [24,25], minimal gaugino medi- 
ation [26], etc. More recently it was pointed out that the same chain may also occur in 
a non-super symmetric context, e.g. Universal Extra Dimensions (UED) [27, 28] and Little 
Higgs theories with T-parity [29]. Therefore, even if the observable SM particles (the quark 
jet and the two leptons) can be uniquely identified, there may still be several competing BSM 
interpretations. Recently there has been a lot of effort on developing various techniques for 
discriminating among different model scenarios [19-22,30-47]. The crux of the problem is 
the fact that the spin of the missing particle A is unknown, and this gives rise to several 
distinct possibilities. Furthermore, the spin of particle A, even if it were known, still does 
not completely fix the spins of the preceding particles C and D. Indeed, since the SM 



^Note that this choice is made only for concreteness of the discussion and does not represent a fundamental 
limitation to our method. AH of our results below can be readily applied in the general case where the visible 
particles are any 3 SM fermions, not necessarily a quark and two leptons. The generalisation of the method to 
the case where the set of visible SM particles includes SM gauge bosons and/or a Higgs boson is straightforward 
and will be presented in a future publication [18]. 
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Table 1: Possible spin configurations of the heavy particles D, C, B and A in the decay chain of 
Fig. |l|. The last column gives one typical SUSY or UED example. In the following we shall use the 
subscript S to label these 6 possibilities. 

particles in Fig. |l| are all spin 1/2 fermions, the particles A, B, C and D must alternate 
between bosons and fermions, but the exact values of their spins are a priori unknown. In 
the spirit of Refs. [21,22], here we shall limit our discussion^ only to particles of spin 1 or 
less, namely we shall consider spin scalars (S), spin 1/2 fermions (F) and spin 1 vector 
particles (V). Table |I] lists the 6 spin configurations for the decay chain of Fig. [l|, which 
were also considered in [21,22]. Five of these six possibilities can be readily accommodated 
in either supersymmetric or UED models. The last column of Table |l| gives some typical 
examples involving the squarks q^ sleptons I and neutralinos Xi in supersymmetry, the KK 
quarks qi, KK leptons £i and KK gauge bosons Zi and 71 in 5D (or 6D) UED [48], and the 
spinless gauge bosons and Zh in 6D UED [49]. The last case in Table I (SFVF) would 
require either a scalar leptoquark or a new gauge boson carrying lepton number. Neverthe- 
less, we include it in our study for completeness and also to connect to the results of [21,22]. 
We should emphasize from the start that we list the supersymmetry and UED examples in 
Table || only as an illustration and in what follows we shall never restrict ourselves to any 
particular model. In particular, we shall not assume any features of the mass spectrum or 
the couplings which might be expected in SUSY or UED. For example, we shall not assume 
a degenerate mass spectrum for the cases which might be expected in UED models, nor shall 
we assume any specific chirality structure of the couplings as predicted in supersymmetry or 
UED. We shall instead keep the spectrum completely arbitrary and also use the most general 
parametrization for the couplings of the heavy partners. Furthermore, we shall not make 
any assumptions about the nature of particle A - it may or may not be the lightest heavy 
partner, and it may or may not be stable. While the dark matter problem mentioned at 
the beginning does provide good theoretical motivation to look for missing energy signals, 
particle A here does not at all have to be the dark matter WIMP, e.g. it may very well decay 
to other heavy particle states, or even directly to SM particles. Consequently, the results 
presented in this paper will be completely general and can be applied to any model of new 
physics which exhibits a decay chain of the type shown in Fig. |l[ 

The main goal of this paper is to assess the possibility of discriminating between the six 

^Our metliod is nevertheless completely general and can be immediately generalised for higher spin particles 
as well. 
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different alternatives in Table ||, using the experimentally observable invariant mass distribu- 
tions of the visible particles (the quark and the two leptons) in Fig. |l|. If such a discrimination 
could be made in a completely model-independent fashion, one could honestly claim a true 
measurement of the spins of the new particles. As a byproduct of our method, we shall also 
obtain an independent measurement of certain combinations of couplings and mixing angles 
of the heavy partners. The invariant mass distributions (of the quark and leptons) are conve- 
nient because they are Lorentz invariant quantities, and are certainly sensitive to the spins of 
the new particles. However, extracting spin, coupling and/or mixing angle information out of 
them is a highly nontrivial task and to the best of our knowledge has not been demonstrated 
up to now in a model-independent setup like ours. The main difficulties can be classified into 
two categories, experimental and theoretical, which we shall now discuss in some detail. 

1.1 Experimental challenges 

This class of problems is related to the ability of the experiment to uniquely identify the 
particles coming from the cascade of Fig. |l|. 

El Jet combinatorics. The events in which the cascade decay of Fig. |l| occurs, will also 
typically contain a number of additional jets. Some of those may come from initial state 
radiation, others may originate from the opposite cascade in the same event, and there 
may also be jets appearing from the decays of heavier particles into particle D. This 
poses a severe combinatorics problem: which one of the many jets in the event is the 
correct one to assign to the D decay in Fig. Some of the existing spin studies in the 
literature simply take for granted that the correct jet can be somehow identified, others 
select the jet by matching to the true quark jet in the event generator output, which is 
of course unobservable. The severity of the jet combinatorics problem is rather model 
dependent and how well it can be dealt with in practice depends on the individual case at 
hand. For example, if the mass splitting between D and C is relatively large, one might 
expect the jet from the D decay to be among the hardest in the event, and this fact can 
be used to improve the purity of the sample. Fortunately, there exists a method (the 
mixed event technique) which should, at least in principle, remove the effect from the 
wrong jet combinations [13]. More recently, the method has been successfully applied 
to measuring SUSY masses at the SPSla study point [50]. A subtraction by a mixed 
event technique is particularly well suited for our purposes, since our method for spin 
measurements only relies on the shapes of the global distributions, and we do not need 
to guess the correct jet on an event by event basis. 

E2 Lepton combinatorics. There is an analogous combinatorics problem related to the 
selection of the two leptons in the cascade of Fig. |l[ First, in general, there may be 
additional isolated leptons in the event, so one might consider requiring two and only 
two leptons per event. However, even then, it is not guaranteed that those two leptons 
are coming from the process in Fig. [l|: for example, each of the two leptons may come 
from a different cascade. Fortunately, there is again a universal method (opposite lepton 
flavor subtraction) which solves both of these lepton combinatorics problems [13]. One 
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forms the linear combination of {e~^e~} + {iJ,~^fi~} — {e~^fj,~} — {/x"'"e~}, in which the 
effects of the uncorrelated leptons in the signal (as well as all SM backgrounds involving 
top quarks, b-jets and W bosons) cancel out^. In what follows we shall be assuming 
that the measured invariant mass distributions have already been properly subtracted 
to take care of the above mentioned jet and lepton combinatorial problems. 

E3 Quark- antiquark jet ambiguity. The cascade shown in Fig. ^ consists of two separate 
processes. In the first one we produce a particle D, which decays to a quark jet and a 
particle C. In the conjugate process, the antiparticle of D is produced and it decays to 
an antiquark jet and the antiparticle of C. Since the two types of jets appear identical 
in the detector^, we cannot distinguish between these two cases, and the observable 
invariant mass distributions are the sum of the individual contributions from these two 
processes. This is a problem since, as we shall see, the sum tends to wash out to some 
extent the spin correlations which may have been originally present. In section ^ we 
shall first present our formulas for the individual quark and antiquark jet distributions, 
but from section ^ onwards we shall always be adding up the quark and antiquark 
contributions together, and we shall use the term "jet" to refer to either a quark or 
an antiquark. For example, when we discuss a "jet-lepton" distribution {ji} we shall 
always imply that it was constructed by adding up the individual quark-lepton and 
antiquark- lepton distributions {q^} + {q^}, so that this quark-antiquark ambiguity does 
not represent a problem. 

E4 Near and far lepton ambiguity. While the charge of the two leptons can be measured 
very well, a priori one does not know which of them is the "near" lepton ^„ (i-e., 
coming from the decay of C) and which is the "far" lepton if (i.e., coming from the 
decay of B). Strictly speaking, once the mass spectrum of A, B, C and D is known, 
one can select a subsample of the original events, in which in and if can be uniquely 
identified. This can be done simply by ordering the two invariant masses mj£+ and 
rriji- as rnjf'^ = max{mj£+, mj^-} and m^°^ = mm{mj£+,mj(-}, and selecting only 
those events for which rnj"^^ happens to be above the observed kinematic endpoint 
of the m}°^ distribution. For that limited sample of events one can unambiguously 
identify i^ and if. However, the price to pay is that the statistics becomes very limited, 
especially if the kinematic endpoints of the rnj'^^ and rii}°^ distributions are close to 
each other. We therefore choose not to apply this trick, and instead we shall consider 
the combined rriji^ and mj£^ distributions for each of the two possible lepton charges. 
This allows us not only to avoid the near-far lepton ambiguity, but also to use the spin 
information contained in the ruj^^ distribution. Previous studies on spin measurements 

^The method is not limited to dilepton events and can also be applied to events with 3 or more leptons. In 
that case one would use all possible dilepton combinations, but include a weight factor for their contribution 
to any given distribution, so that the total weight of any given event, summed over all dilepton combinations, 
is 1. 

*If g is a heavy flavor, the distinction can be made (statistically). To be conservative, we ignore this 
possibility in order to demonstrate that our method works even in the worst case scenario of jet ambiguity. 
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have concentrated on the spin correlations between the jet and the near lepton, for 
which relatively simple and compact analytical expressions can be derived. The jet-far 
lepton contribution was regarded to a large extent as an annoying background which 
tends to wash out the jet-near lepton correlations. Our approach is very different: we 
actually treat both rriji^ and rriji^. distributions on the same footing. Since we have 
derived the most general expressions for both rriji^ and mj£^, in our method we are in 
effect able to fit separately to each one, and we do not even need to make the 
discrimination on an event by event basis. In this sense our method is using all of the 
available information about spins which is present in the data. 

Additionally, there are the usual complications on the experimental side, such as SM back- 
grounds, detector acceptance and resolution, triggering etc. All of these factors should be 
taken into account when trying to decide how well our method will work in any particular 
case. But the main advantage of our method is that it is completely general, and can always 
be applied, even in the extremely complex environment of a hadron collider experiment. 

1.2 Theoretical issues 

Even if none of the experimental issues E1-E4 discussed above ever existed, e.g. we had a 
perfect detector, and we could somehow identify on an event by event basis with absolute 
certainty which particular jet and two leptons came from the cascade in Fig. ^, and further- 
more, we could discriminate q from q as well as £n from even in that idealized case, there 
would still have been a long way to go towards a clean spin measurement, i.e. a discrim- 
ination between the 6 cases of Table |l|. The problem is that the measured invariant mass 
distributions depend on all of the following 4 factors: 

Tl Mass spectrum. It is well known that the shapes of the observed invariant mass distri- 
butions in general depend on the heavy partner spectrum. In fact this has been used in 
the past to make mass measurements of the heavy partner masses, especially in the case 
when one of the heavy particles in the chain is off-shell [51,52]. Mass measurements 
are therefore a useful (but not necessary - see below) first step towards determining 
the spins. For simplicity, throughout this paper we assume that all masses m^, rns, 
mc and mo have already been determined from kinematic endpoints. This assumption 
is common with all previous spin studies. It appears rather feasible, since the mass 
measurements only require the extraction of the kinematic endpoints, which are sharp 
features in the invariant mass distributions, and those are likely to be seen in the data 
much earlier than the actual shape of the distributions. However, we should emphasize 
that our assumption about the known mass spectrum was made only for simplicity, and 
to keep the discussion focused on the more challenging measurements like the spins, 
couplings and mixing angles. Our method in fact does not require any prior knowledge 
of the mass spectrum. When the mass spectrum is a priori unknown, the fits described 
in Sec. ^ would actually pick up the correct values of the masses, in addition to the spin 
and coupling measurements. 
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T2 P article- antiparticle ambiguity (D/D). This problem is related to the experimental 
issue E3 from the previous subsection. Since we do not know if the jet was initiated by 
a quark or an antiquark, we also do not know whether the heavy particle cascade was 
initiated by a particle D or its antiparticle D. At a pp collider such as the Tevatron, 
the symmetry of the initial state implies that the fraction f of D particles produced in 
the data should be equal to the fraction / of antiparticles D. Unfortunately, at a pp 
collider like the LHC, the initial state is not symmetric, so one may expect an excess of 
particles over antiparticles: / > /, but the precise value of this excess A/ = / — / is a 
priori unknown. Therefore at the LHC / is in principle an unknown parameter, which 
significantly affects the observable and invariant mass distributions. Most 

previous studies of spin measurements have fixed / to the value for the corresponding 
study point [19,20]. However, in the absence of an independent measurement of /, this 
is unjustified. The influence of / on the spin extraction was considered in [38,40], where 
/ was left as a floating parameter and consequently the extraction of the spins became 
much more difficult. In what follows we shall follow a similar approach, namely, we 
shall not make any assumptions about the value of / when we discuss measurements at 
the LHC and we shall instead treat / as a free input parameter. Only in Sec. [4.2| , where 
we apply our method to the Tevatron, we shall take f = f ■ Naturally, / is trivially 
related to / as 

/ + / = 1. (1.3) 

T3 Chirality of the fermion couplings. Note that the three SM particles in Fig. ^ are all 
fermions, whose couplings to the heavy partners at each vertex are a priori unknown. 
The observed invariant mass distributions depend on the chirality of those couplings, 
and this presents a formidable challenge in measuring the spins. The problem is that 
any given set of measured invariant mass distributions could in principle be explained 
by one spin configuration with a certain choice of chiralities, or a different spin config- 
uration with a different choice of chiralities for the fermion couplings. To the best of 
our knowledge, none of the existing spin studies have accounted for this ambiguity in 
a consistent and fully model-independent way. Our main objective in this paper is to 
devise a method for spin measurements which makes no assumptions about the chirality 
of the couplings at each vertex in Fig. |^. Correspondingly, we shall keep those couplings 
completely arbitrary, and parameterize them in the most general way in terms of inde- 
pendent chirality coefficients at each vertex. For example, in the case of an interaction 
between a heavy spin 1/2 fermion F, a heavy scalar $ and a SM fermion / we take the 
interaction Lagrangian to be 

C{F,f,^) = ^F{9LPL + 9RPR)-^f^ + h.c. (1.4) 

where gi and gn are arbitrary (and in general complex) coefficients. In general, there 
are three different sets of {gLidR]: one at each vertex of Fig. |l|. We shall denote them as 
{cliCr}., and {aL,a^}, as shown in Fig. ||. Similarly, in case of an interaction 

between a heavy spin 1/2 fermion F, a heavy vector boson and a SM fermion / we 
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use the interaction Lagrangian 



CiFJ,A^) = ^Fl^{gLPL+gRPR)^fA^ + h.c. (1.5) 

where just Uke before the coefficients {gL,gR} = {cl,cr}, {6l,6_r} or {aL,aR}, de- 
pending on the vertex. In what follows we present our results in terms of these most 
general coefficients {cl,cr}, {bL,bji} and {gljCir}. According to our convention, the 
couplings {cl,cr} are always associated with the D-C-q vertex, the couplings 
are always associated with the C-B-ln vertex, and the couplings {0^,0/?} are always 
associated with the B-A-if vertex. We shall not be specifying explicitly whether a 
given pair such as {aL,aR} parameterizes the interaction (|1.4|) or the interaction ( |1.5| ), 
since that should be clear from the context. 

We shall see below that the shapes of the invariant mass distributions only depend 
on the relative chirality of each vertex, therefore it is convenient to unit normalize the 
couplings as 

IolI' + |ai?P = 1 , (1.6) 
+ = 1 , (1.7) 

\cl\^ + \cr\^ = 1 , (1.8) 

In that case, the relative chirality at each vertex is parameterized in terms of a single 
parameter, which can be taken as an angle: 

I^-rI \^r\ |c_r| 

tanipa = -. r, tanipb = -rr—r , tanipc = -. — r. (1.9) 

PlI \cl\ 



2 J 



By convention, we shall take all three of these angles to be defined in the range [0, 
(as opposed to [vr, ^]). The angles ipa, and ipc encode all of the relevant^ model 
dependence, e.g. the nature of the interaction and the mixing angles of the heavy 
partner mass eigenstates. It is worth emphasizing that we consider the couplings gL 
and gR in eqs. (|1.4 |1.5D to be the couplings in the mass eigenstate basis for the heavy 
partners. Therefore, whenever there is mixing among the heavy partner states, our 
couplings gL and gR are in general matrices which are related to the couplings g^^^ and 
in the interaction eigenstate basis through rotations by the corresponding mixing 
angles 

9L,R = Uf^ Ub , (1.10) 



^At this point it may be useful to do a quick count of the relevant degrees of freedom. For example, consider 
the B-A-lf vertex parameterized by {aL,aR}- Since sl = |a_L|e'^^ and an = |ai{|e'*" are in general complex 
parameters, originally there are four degrees of freedom {\aL\, \a.R\, 4>l and (jya) parameterizing each of the 
SM fermio n in teractions (1.4 Ls]) . One combination of \aL\ and \aii\ is eliminated through the normalisation 



condition (L6), while (L9) simply parameterizes the other combination of \aL\ and |a_R| in terms of ipa- The 
remaining two degrees of freedom, the phases (j>L and <^h, remain arbitrary and cannot be measured from the 
invariant mass distributions that we are considering here. Instead, they will have to be measured by some 
other means. 
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where the matrix Up (Ub) diagonahses the mass matrix of the corresponding heavy 
fermion (boson). Due to this mixing, in general we do not expect our couphngs ql and 
qr to be purely chiral, even in models where one starts with purely chiral couplings 
g^^"^ and g^^^ in the interaction eigenstate basis. The effect of heavy fermion mixing Uf 
in a specific UED model was previously considered in [41], and here we generalise the 
discussion to the case of arbitrary heavy fermion mixing Uf, arbitrary heavy boson 
mixing Ub, and arbitrary couplings g^^^ and g^^\ Clearly, there is an enormous num- 
ber of model-dependent parameters contained in Uf, Ub, g^^ and g^^\ and it will be 
rather hopeless to try to measure them all at once. One of the main results of this paper 
will be to identify which particular combinations of these coupling and mixing angle 
parameters can be experimentally measured from the invariant mass distributions of 
the three SM fermions (in our case, q, in and if), and to propose the actual method for 
measuring them. We shall find that there are three such combinations, which we shall 
call a, P and 7 (for details, see Sees. ^ and 5^). Each one of them is potentially exper- 



imentally accessible, and represents some combination of couplings and mixing angles 
as illustrated in eq. ( |1.10| ). It is in this sense that our method yields a measurement of 
the couplings and mixing angles of the heavy partners, as advertised in the abstract. 

T4 Spins. Finally, the invariant mass distributions also contain information about the spins 
of the heavy particles along the decay chain. For example, pure phase space predicts 
flat (in m^) invariant mass distributions for SM particle pairs originating from adjacent 
vertices in the decay chain. Deviation from this pure phase space prediction implies 
some kind of spin correlations [19] , but what type? Conversely, observing distributions 
which are consistent with the pure phase space prediction does not necessarily mean 
that all particles involved in the decay are scalars - spin correlations may have been 
present for the individual subprocesses (to be defined below) but may have been washed 
out when added up to form the experimentally observable distributions. Below we shall 
encounter examples of both of these situations. 

The general approach in previous spin studies has been to compare the data from a given 
study point within one specific model to the corresponding data obtained from another model 
alternative with different choice of spins for the heavy partners. A common flaw in all such 
studies was that three of the four relevant factors, namely Tl, T2 and T3, were fixed to 
be identical in the two models, so that any remaining difference can be interpreted as a 
manifestation of spins (the factor T4 above). However, this is not the correct approach when 
it comes to actual pure measurements of spins in a model-independent fashion. Since the 
chirality parameters (pa, '■Pb and c/^c and the particle-antiparticle ratio / are not independently 
measured prior to the attempted spin determination, they need not have the same values for 
each of the different spin configurations under study (in our case, the 6 ones listed in Table |^) 
and should be allowed to float. Therefore, the proper question to ask instead is: 

Given the data, which (and how many) spin configuration gives a good fit to it 
for some choice of the chirality parameters (pa, 'Pb and ipc, and for some choice of 
the particle-antiparticle ratio /? 
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The main result of this paper is that we provide the tools needed to address this question 
in a completely model-independent way, namely in order to determine whether a given spin 
configuration "S" is consistent with the data or not, we do not need to specify the values of 
/ and /, nor do we need to specify the chirality of the couplings (/9a, ^ph and ipc- In other 
words, we have divided the question posed above into two parts: for a given mass spectrum 
(i.e. factor Tl is known), 

• Ql: What is the spin, i.e. what is factor T4? 

• Q2: What are the particle-antiparticle fractions / and / (item T2 above) and what are 
the couplings and mixing angles (item T3 above)? 

Our method allows us to provide an independent answer to the spin question Ql regardless of 
the answer to the follow up question Q2. In this sense we are able to make a pure measurement 
of spin in a model-independent way. Of course, as we shall see below, the actual answer to 
the question Ql may not be unique, and sometimes there are cases where more than one 



particular spin configuration may fit the data. In fact in Sec ^Jj we shall show that the 
model pairs {FSFS, FSFV} as well as {FVFS, FVFV} are quite often indistinguishable. 

Since we have decoupled the spin issue T4 from the /-/ issue T2, our method is not 
limited to pp colliders such as LHC, and is equally applicable to the Tevatron. In contrast, 
the lepton charge asymmetry proposed by Barr [19] is greatly affected by the value of /, 
for example it is predicted to be identically zero at the Tevatron and has no discriminating 
power there with regards to spins. In this sense our method provides a pure measurement of 
the spins and the spins alone. What is more, in the process of answering the spin question, 
we also get a measurement of some combination of the couplings and / and /. In this sense 
our method is also the first and most general attempt to measure mixing angles of heavy 
partners (e.g. superpartners) at the LHC. 

The paper is organised as follows. In Section |^ we describe the main idea of our method 
and derive the main building blocks for the spin measurement. In particular, we give ex- 
act analytical expressions for all relevant invariant mass distributions (including {g^j } and 
{9^/ }) in most general case of arbitrary couplings, arbitrary / and /, and arbitrary mass 
spectrum, for each of the six cases from Table ||. Our results in Sec. § generalize those of 
Refs. [17,20-22]. In Sec. ^ we reorganise our results from Sec. § to form the experimentally 
observable invariant mass distributions and We also derive the exact 

combinations of couplings and mixing angles which are being measured as a byproduct of the 
spin measurement^. Section ^ begins by summarising the key analytical results from the pre- 



vious two sections, and outlines our method for spin and coupling measurements. In Sec. 4.1 
we prove analytically the degeneracy of the {FSFS, FSFV} and {FVFS, FVFV} model pairs 
- we derive the relation between the couplings and mixing angles within each pair of models 
which would result in identical observable invariant mass distributions for those model pairs. 



^Readers who are only interested in the practical applications of our results, and would prefer to skip these 
mathematical derivations, are invited to jump directly to Sees. ^ and |5|, which are self-contained and can be 



read independently from the more technical sections ^ and | 
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In Sec. 4.2 we specify our results to the case of pp colliders such as the Tevatron and show 
that our spin analysis method can be just as successful there. Finally, in Sec. ^ we provide 
an illustration of an actual idealised measurement, using a mass spectrum and couplings 
as for the SPSla study point in supersymmetry. Assuming that the data comes from each 
one of the 6 models from Table |^ in turn, we then demonstrate how well the remaining 5 
possibilities can be ruled in or out. This results in a total of 36 different case studies, the 
results of which are presented and analysed in that section. In Sec. ^ we summarize our main 
conclusions, and discuss the pros and cons of our method in comparison to other proposals 
for spin measurements in the literature. 

2. General expressions for the invariant mass distributions 
2.1 Preliminaries 

The basic idea behind our method is the following. For any given spin configuration S, we 
write the invariant mass distribution of a pair of SM particles from Fig. || as 

2 2 

(S2) =Y.Y.Kfhf^^Mc)r^UK-^-^y^^) ^ (2-1) 

where the index p denotes one of the five possible SM particle pairs: p = {jln ' i^n ' i^/ ' i^/ ' }; 



rhp is the unit-normalised invariant mass 



< mp < 1 , (2.2) 



^max 



i.e. the invariant mass nip scaled by the value of the corresponding kinematic endpoint 
which has already been measured from the corresponding nip distribution. The mass ratios 



X, y and z were already defined in (L2), while {U} is a pair of indices denoting one out of 



four possible classes of subprocesses Pjj which will be discussed in detail below in Sec. 



The coefficients K^^j and the functions will be explicitly defined later in Sec. 2.2 



The general expression (2.1) corroborates our discussion in Sec. 1.2 ~ we see that the 

invariant mass distributions indeed depend simultaneously on all of the four factors (T1-T4) 

discussed earlier. However, notice that the coefficients Kjj in the expansion (^]^) only depend 

on the particle/antiparticle fraction / and the chiralities ipa-, and (pc, i.e. factors T2 and 

T3. On the other hand, the functions J-g!\j{rh?;x,y,z) only depend on the mass spectrum 

(factor Tl) and the spin (factor T4). Once the spectrum is measured and the mass ratios 

(p) 

X, y and z become known, the functions J'g.jj only depend on rh and provide a unique basis 
which can be fitted to the data for each of the measured distributions {p} . Since the functions 
do not depend on the model dependent parameters /, ifa-, ^fb and tpc, this fit can be done 
in a completely model-independent way, without any prior knowledge about the nature of 
the particles A, B, C and D, the nature of their couplings, or the size of their mixing angles. 
For each of the 6 possible spin configurations S, this fit may or may not yield a good match: 
then, those spin configurations which give a bad fit to the data will be ruled out. Conversely, 
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the spin configurations which give a good fit will be ruled in, and furthermore, the values of 
the fitted coefficients K will represent a measurement of the couplings and mixing angles of 
the heavy partners. 

2.2 Classification of helicity combinations 

Table |2| lists all possible helicity combinations (32 altogether) contributing to the process 
of Fig. |l|. The 8 combinations shown in blue have been previously considered in [20-22]. 
The remaining 24 combinations shown in red are being considered here for the first time. 
We find it convenient to classify all possibilities into four categories Pij,I,J = 1,2, where 
each category gives rise to the same functional dependence for the three invariant mass 
distributions of interest: {j^J}, {j-^/} and {^^^J}- We name these four categories as follows: 

• Processes of type Pu. These include all cases where the physical helicities of the 
(anti)quark jet and near lepton are the same, while the physical helicities of the two 
leptons are opposite. The four processes of type 1 in the nomenclature of Refs. [20-22] 
fall into this set. In addition in this group we find four new combinations involving 
right-handed quarks. 

• Processes of type P2i- These include all cases where the physical helicities of the 
(anti)quark jet and near lepton as well as the physical helicities of the two leptons 
are opposite. The four processes of type 2 in the nomenclature of Refs. [20-22] fall 
into this set. Again, there are four new cases involving right-handed quarks. Note that 
the processes of type P21 are simply obtained from those of type Pu by interchanging 
q ^ q while keeping the chirality labels fixed. 

• Processes of type P12. Here the physical helicities of the (anti)quark jet and near lepton 
as well as the physical helicities of the two leptons are the same. These processes are 
obtained from those of Pu by changing the chirality label of the far lepton: L <^ R for 

• Processes of type P22- Here the physical helicities of the (anti)quark jet and near 
lepton are opposite, while the physical helicities of the two leptons are the same. These 
processes can be obtained from P12 by interchanging q <^ q, or alternatively, from P21 
by changing the chirality label of the far lepton: L <^ R for 

All processes falling into the last two categories are new, and more importantly, as we shall see 
below, they give a qualitatively new functional dependence of the dilepton and ji j invariant 
mass distributions which was not exhibited in the previous studies [20-22] . 

It is worth noting that in the case of a heavy fermion (F) , there is a distinction between the 
Dirac and Majorana case. For a Dirac fermion, half of the processes within each category Pjj 
of Table |2| are absent, since the adjacent SM fermions must be a particle and an antiparticle. 
For a Majorana fermion, there is no such restriction, and all processes exhibited in Table |2| 
are in principle allowed. 
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Processes Pu 


Processes P12 




/IcLpl^LplaLp 


f\cL\'MaR\'' 


{qlJI^^r} 


f\cL\'\bRnaR\' 


flcLl'MaRl^ 


fMbR^aLl' 


/|cLp|6i?P|ai|2 


f\cR\'\bR\^\aii\^ 


r — /?+ /? — 1 

/|cflP|6R|2|aR|2 


{qr^rJU 

flcRl^MaLl' 


r — /?+ /) — 1 

{qrJbJl} 
f\cR\^\bR\''\aL\'' 


f — n — /?+ ~i 

{(iRAL^q} 

]\cR?\bL?\aL? 


{qR^lJi} 

f\cR\^\bL\W 


{^R.^L^^r} 
f\cR?\bL?\aR\^ 


Ur^H^^r} 
flcRl^MaRl'' 


f\cL\^\bL\'\aL\' 




]\cL?\bL?\aR\^ 


f\cL\'\bL?\aR\^ 


f\cL\'\bR\'\aR\' 


f\cL\^\bR\^\aR\^ 


f\cL\'\bR\^\aL\' 


{Ql, 1^,11} 
f\cL\^\bR\'\aL\' 


{QR^^w^r} 
f\cR\^\bR\^\aR\^ 


{QR^^w^'r} 
f\cR\^\bR\^\an\^ 


{qRj-RJt} 

f\cR\''\bR\^\aL\^ 


UrJIiJI} 

f\cR\^\bR\''\aL\'' 


{QR^l,^!} 
flcRl'^MaLl^ 


{QRJt^^l} 
/|ci?j2|6ip|aip 


{qr^l^^r} 
f\cR\^Man\' 


/|cifj2|6z,p|aR|2 


Processes P21 


Processes P22 



Table 2: Classification of all possible hclicity combinations contributing to the process of Fig. |l|. 
The combinations shown in blue have been previously considered in [20-22] . The combinations shown 
in red are being considered here for the first time. Under each helicity combination, we also show the 
associated prefactor contributing to Kj'^J in eq. 



2.3 Invariant mass distributions 

In principle, there are 9 invariant mass distributions that we can form: 

2 2 



dN 



qtn / s 
dN \ 



dmi , ] 

qin / s 

dN 



dm? ^ 



f / s 



dN 



dN 



drh?n 



1 



1=1 j=i 
2 2 



1=1 j=i 
2 2 



1=1 j=i 
2 2 



^ E E ' V'a, V'fe, ^c)J^s-!/J {m'^gf,± ; a;, 2) , 



1=1 J=i 
2 2 



(2.3) 
(2.4) 

(2.5) 

(2.6) 

(2.7) 



1=1 j=i 



where the factor of i on the right hand side was introduced for future convenience. Note 
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that it is the same set of functions J'g.fj which enter both the {qin} and {qin} distributions 

(m^ X, y, z) = (^'; X, y, z) = J'iljiim^; x, y, z) , (2.8) 

and similarly, it is the same set of functions J^g./j which enter the {q(f} and {q£f} distribu- 
tions: 

X, y, z) = rf/]{w?;x, y, z) ^ T^Jjji^^' ^' ^) ■ (2-9) 

(p) 

In the following two subsections we shall separately define and discuss the functions jj 
and the coefficients Kj^J appearing in the general expressions (2.3-2.7). 

2.3.1 The functions J^j?]j 



Eqs. (2.3-2?7|) show that all invariant mass distributions can be written in terms of three 
sets of basis functions: J-'g.f^J {m'^ ; x , y , z) , t'^J.^/J {m^ ■, x , y , z) and T^^^] j{'m? ■, x , y , z) . We shall 
define the basis functions to be unit normalized: 

I J's.u{m'^;x,y,z)dm'^ = 1 , (2.10) 
Jo 



J'[^^/hm^-,x,y,z)drn^ = 1 , (2.11) 



oo 



x, y, z) dm^ = 1 . (2.12) 







With this normalisation, all basis functions J^^-lj {^"^'^x, y, z), J^^^jj{m'^; x, y, z) and ^g.^/J {rn?-^ x, y, z) 
are defined in Appendix 

(p) 

A few comments regarding the J'g jj functions are in order. Recall that half of the pro- 
cesses belonging to category Pu and P21 (in the classification of Sec. 2.2) have been previously 
considered in [20-22], so that the functions J'g'li and ^g'li in principle already appear there. 

We find agreement with [20-22] for the case of J'g'li and J'g'li^ ^^'-^ supplement those 

(p) (p) 

results with the remaining two types of functions J^gi2 •^322- shall now comment 

ip) 

individually on each type p of basis functions jj- 

Table ^ in Appendix |A| shows that the J'gj'j functions are pairwise equal: 

J^^lll\m^;x,y,z) = J^f{^ {rn'-x^y, z) , (2.13) 

J^f^^{m'';x,y,z) = J^f^:^ {rn^-x,y, z) . (2.14) 

These relations are easy to understand: processes P12 differ from processes Pu only by the 
chirality label of the far lepton If. However, the j£„ distribution does not know about the far 
lepton, therefore the J'g./j function should be the same for both J = 1 and J = 2. Table |6| 
has essentially already appeared in [21] (see Tables 10 and 11) and we reproduce it here just 
for completeness. 



- 15 - 



On the other hand, Table ^ of Appendix ^ contains some new results for the ^^g^jj 
functions. In this case there are still only two independent functions, but the functional 
relationship is different from ( |2. 13 , 2.14): 



r^^^^^{fn?;x,y,z) = J^^g^^^{m^;x,y,z) , (2.15) 

•^s!i2("^^; ^) = ^s%^^'^'^ ^) • (^-i^) 

Again, the reason behind these relations is easy to understand intuitively. Processes Pij are 
related to processes P2J by simply interchanging q ^ q, which, of course, does not affect the 



two leptons which are further down the cascade decay chain. Because of ( 2.15 ), Refs. [20-22] 
found identical results for J-g^ii and J's2i (corresponding to processes of type 1 and 2 in 
their notation), but missed the functions ^'^^12 and ^'^'22- This was a direct consequence of 
the underlying model dependence, and in particular factor T3: the studies [20-22] assumed 
very specific fixed values of the chirality coefficients (namely, cl = 1, cr = 0, bi = 0, 5^ = 1, 
a/, = 0, Or = 1 for the supersymmetry example and cl = 1, cr = 0, bi = 1, bji = 0, ai = 1, 
an = for the UED example) and therefore their results, while correct, are only valid within 
this limited model-dependent context. In contrast, deriving the complete set of functions 
•^S% ^'^^ possible sets of processes Pjj allows us to address the spin question Ql raised 
in the Introduction in a completely model-independent fashion. 

Similar remarks hold for the J'g.^j functions in Appendix ^ Here again the functions 

•^s-ri ^^'^ -^3-21 ^S^ce'^ with the results of [21], while the functions ^^.1^2 ^^'^ 
new. However, whether (and what type of) relations exist between the four functions J'g./j 
varies from case to case (i.e. the value of the spin configuration index S). In the three cases 
(SFSF, FSFS and FSFV) where there is an intermediate heavy scalar between the emitted 
jet and far lepton, the J^g./j set is again reduced to only two independent functions, however, 
the exact functional relations are also S'-dependent: for S = 1 (SFSF) we find 

J^aii^ (m^; X, y, z) = T[ili> {m'; x, y, z) , (2.17) 

^g/^ (m2; X, y, z) = (jr?; x, y, z) , (2.18) 

while for 5 = 2 (FSFS) and 5 = 3 (FSFV) we find 

^s,ii(.^^'^ X, y, z) = J^^ll{\m^;x, y, z) for 5 = 2, 3 , (2.19) 

J^^-J'j^^(m2;x,y,z) = J^^^'2^^(m2;x,y,z) for 5 = 2,3 . (2.20) 

In the remaining 3 cases S = 4,5,6 (i.e. FVFS, FVFV and SFVF) we find that all four 
functions J^^g\^] are independent. 



^The only discrepancy we found was in the constant coefficient in front of the \ny and Inm'^ terms in the 
J-'^-iP function: in eq. (B.9) of Ref. [21] it is listed as —(2 + Ay) while we find —(1 +4jy)2. Since our results 
agree with the numerical results of Figs. 5a and 5b in [21], we believe that eq. (B.9) in [21] has a typo. 
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2.3.2 The coefficients K 



ip) 
IJ 



Having defined the complete sets of functions J'^'-ij entering the general expressions (|2^ 

|2.7| ), it now remains to define the coefficients Kj^j{f]ipa,(pb,^c) entering those formulas. 
Notice that these coefficients do not carry a spin index S, i.e. they are independent of the 
assumed spin configuration. Therefore we only need to define them for each fermion pair 

Using the factors from Table for the coefficients belonging to processes Pn we readily 
obtain 



K 
K 
K 
K 



11 

11 

11 

(get) 
11 



f\cL\^\bL\^\aL\^ + f\cR\^\bR\^\aR\\ 



7) 



^11 
-'^11 



L\ |«L| 

f\cL\'\bR\'\aR\' + f\cR\'\bL\'\aL\' 
f\cL\'\bL\'\aL\' + f\cR\^\bR\'\aR\^ 



(2.21) 
(2.22) 
(2.23) 
(2.24) 



The corresponding coefficients for processes P12 can be now simply obtained from (2.21-2.24) 
by the substitution ai <-> qr: 



K 



K 



K 



K 



12 

12 
12 
12 



K 
K 



«) 
12 

(9-4) 
12 

{qr,) 



12 



-"■12 



f\cL?\bL?\aR\'' + f\cR\''\bR\^\aL\\ 

f\cL?\bR\^\aL? + f\cR\^\bL\>R\\ 

f\cL?\bR\^\aL? + f\cR\^\bL?\aR\\ 
f\cL?\bLWR\^ + f\cR\^\bR\^\aL\\ 



(2.25) 
(2.26) 
(2.27) 
(2.28) 



Next, replacing f f and q ^ q in (2.21- 23^ ) gives the corresponding coefficients for 
processes P21: 



K. 



{q^n) 

21 



K. 



K. 



21 

t) 

21 
21 



^21 

-"^21 
-"■21 

-'^21 



f\cL?\bL?\aL? + f\cR\^\bR\''\aR\\ 

f\cL?\bRWR\'' + f\cR\^\bLWL? ^ 

f\cL?\bR\^\aR\^ + MHWl? ■ 

f\cL?\bL\>L? + f\cR\^\bR\^\aR\\ 



Finally, replacing a/j in (2.29- 23^ ) yields the coefficients for processes P22: 



K. 



22 

(qi^n) 

22 



K. 



(9-4) 



22 



(qip 



K. 



K. 



22 

{qi^t) 
22 



^22 
-"■22 

-^^22 



f\cL?\bL?\aR\^ + f\cR\^\bR\\L? ^ 
f\cL\^\bRWL? + f\cR\''\bLWR\^ ^ 
f\cL?\bR\^\aL\^ + f\cR\^\bLWR\^ . 
f\cL?\bL\>R? + f\cR\^\bR\^\aL\^ . 



(2.29) 
(2.30) 
(2.31) 
(2.32) 

(2.33) 
(2.34) 
(2.35) 
(2.36) 
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The coefficients Kf^ for the dilepton distributions can be expressed in various ways, for 



example in terms of the coefficients involving the near lepton 



T^W _ r^(ql-) r^(ql-) r^{ql+) r^(ql+) . 
^IJ —^IJ ~^^IJ ~^^IJ '^■^IJ ' 



in terms of the coefficients involving the far lepton if. 



K) 



(g£7) (q£l) (qij) (qij) 



in terms of the coefficients involving the positively charged lepton 



K 



or finally, in terms of the coefficients involving the negatively charged lepton 



(2.37) 



(2.38) 



(2.39) 



(2.40) 



All of the definitions ( 2.37 - 2740| ) are equivalent because of the relations ( 2.21 - 2l3^ ) existing 
between the various coefficients. Notice the normalisation condition 



1=1 ,7=1 

With the definitions (2.21- 2l3^ ) and the conventions ( |2.1C - p.l^ ) and ( |1.6Hl.8| ), our distri- 



butions (2.3-2.7) are normalised as follows: 



dN 



drh? 



qtn / S 



dN 



dN 
dmj, 



drni ± 



dm„ 




drh? i,± 



drh?_.± 



f_ 

2 ' 



2 ' 



1 . 



(2.42) 

(2.43) 
(2.44) 



It is now clear how the factor of ^ in eqs. (2^2.7) is related to the normalisation: the 
dilepton distribution (^]^), which is experimentally observable, is unit normalised, as seen by 
eq. (|2.44|) . On the other hand, eqs. (|2.42|) and ( |2.43|) show that the individual [qin], {qin], 
{qif} and {qif} distributions are not unit normalised. However, this is not a problem, 
since those distributions cannot be separately observed. In fact, as we shall see in the next 
section, the normalisation ( p.42 , 2743| ) is precisely what is needed in order to unit normalise 
the observable invariant mass distributions for {ji^} and 
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3. Observable distributions in a {q^i'^ji^} chain 
3.1 Invariant mass formulas in the {J^g'lj} basis 

If we could identify the nature of the jet {q versus q) on an event by event basis, we could 
use directly the distributions (|2.3|-2.7) derived in the previous section. As mentioned in the 



Introduction, there may be cases where this is possible, e.g. if g is a 6-quark, or alternatively, 
if it is a lepton so that the decay chain of Fig. |l] represents a trilepton signature. Here, 
however, we shall make the conservative assumption, which also happens to be true in many 
models, that g is a light flavor quark, so that the experimental distinction between a q and q 
cannot be made. In that case, we have to add the corresponding distributions involving a q 
and a q: 

f dN \ f dN \ f dN \ 



+ 



d7v?.„+ / \ drh?„+ / \ dim?_„+ j 

J^n / S ^ l^n / S ^ l^n / S 

2 2 

= li2i2^if^(f^'P-'^b,^c)4fj\^'e^-^^^y^^) ' (3-1) 

7=1 J=l 

dN \ I dN \ I dN 

+ 



J'-f / S \ ^ f / s \ ^ < 
2 2 , . 



^ \iZil^T/\f^^a,^,,^c)J'^i'/]{rhl±-,x,y,z) . (3.2) 
/=i j=i ^ 

Since the functions do not depend on the q-q ambiguity (factor E3), the new set 

of coefficients Kfj" and -f^/j can be simply related to those already introduced in the 
previous section: 

K^'j"\f, ^a, ^b, Vc) = K\f\f, Va, ^b, ^c) + Kf;-\f, Vb, ^c) , (3.3) 
Ku^\f, ^a, V>b, V>c) = K'fj^\f, ipa, ^Pb, ^c) + K^ij'\f, ^Pa, V>b, V>c) ■ (3.4) 



Substituting the definitions ( |2.21 - ^^ ) into ( |3.3| ) and (3.4), we find that the Kj'^j^ coefficients 



can be expressed in terms of the particle-antiparticle fraction / and the relative chiralities 

ipa, ^b and ipc as follows 







V'b, 




= {f\cL 


^ + f\cR 




aL 


' + (/>L 


' + f\cR 


')\bR\' 


aR\^ , 


(3.5) 


-"■12 


\f, ^a, 


V'b, 


"Pc) 


= {f\cL 


^ + I\CR 




aR 


' + {f\cL 


' + f\cR 


')\bR\' 




(3.6) 


-"■21 


\f, ^a, 


V'b, 


"Pc) 


= {f\cL 


^ + f\cR 




aL 


' + {f\cL 


' + f\cR 


')\bR\' 


aiJp , 


(3.7) 


^22 


\f,^a, 


V^b, 


'Pc) 


= {f\cL 


^ + f\cR 




aR 


' + {f\cL 


' + f\cR 


')\bR\' 




(3.8) 
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Figure 2: A contowc plot of cos(fc as a function of cos (pc and f. 
The remaining Kfj coefficients can be related to these as 



^11 


— -"-11 


— -"-21 


— -"-21 




— ^12 


_ T^U^t) 
— -"-22 


— -"-22 


^21 


— ^21 
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(3.9) 
(3.10) 
(3.11) 
(3.12) 



It is important to notice that while the coefficients Kj''j\f,ipa,(pb,^c) defined in ( 3l5| - 3.1 j ) 
depend on all four variables /, ipa, (fb and ipc, the dependence on / and ipc only appears 
through the combinations /|clP + /|c/jp = f cos"^ ipc + fsin^ipc and /|clP + /|c_rP = 
/ cos^ (^c + /sin^ ipc- We shall therefore find it convenient to introduce an alternative chirality 
parameter (fc defined by the relations: 



so that 



cos^ (^c = / cos^ <fc + f sin^ (fc , 



sin^ Cpc = f cos^ <^c + / sin^ <^c , 



cos 2(pc = if - f) cos 2ipc 



(3.13) 
(3.14) 

(3.15) 



The relationship between the newly introduced parameter ipc and the original parameters / 
and ipc is pictorially illustrated in Fig. ||. 
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In terms of the new parameter (^c, the defining equations (^-3.8) for the K^jP (/, v^ai 'Ph, H^c 
coefficients simply become 
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(3.19) 



and the remaining relations (|3.9H3.12| ) are unchanged. 

Using the relations (3.1(:-0^), and the normalisation conditions (1.3) and (1.6-1.8), it 



is easy to check that the Kfj coefficients obey the following normalisation conditions 

2 2 



1=1 j=i 

2 2 
7=1 J=l 



IJ 



IJ 



1 



1 . 



(3.20) 
(3.21) 



Given the unit normalisation ( |2.1C| - |2?T^ ) of our basis functions J'g^lj, eqs. ( 3.20| ) and ( 3.21| ) 
readily imply that the {jin} ^^^^ i-'^/ J" distributions (|3.lD and (3^) are automatically half- 
unit normalised^ 



poo 


f dN 


I 




n 


^ dN 







(3.22) 
(3.23) 



/ / s 



The last step in deriving the experimentally observable invariant mass distributions is to 
recall that the near and far lepton (^„ and £/) cannot be distinguished on an event by event 
basis, therefore we need to form the distributions which are based on definite lepton charge: 



dN \ 



dN 



dm?., j 
]^ / s 



dN \ 


/I 


( dN 






dN \ 




( dN 




dm^„ 



(3.24) 



(3.25) 



*This can also be seen directly from the definitions ( |3.f[ ) and (3.2) of the j£ distributions and making use 
of eqs. ( p!42| ), ( plsj ) and 
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When combining the jet-near lepton and the jet-far lepton distributions in eqs. ( p.24| , ^35| ), one 
has to be careful since until now each individual distribution was written in terms of its own 
unit-normalised invariant mass variable rhji^ and rhjij, . In general, these two variables will be 
different, since the kinematic endpoints mlj^^ and rn^^^ , to which they are normalised, will 
not coincide. Once this problem is identified, it can be handled in various ways, for example, 
by writing out the sums ( p .24 , ^^251 ) in terms of the actual (i.e., not unit-normalised) invariant 



masses. In this paper, we prefer to keep the fh notation, and write all of our distributions in 
terms of unit-normalised invariant mass variables. To this end, we normalise any jet-lepton 
invariant mass mj£ to the endpoint 

= max{m,"^f , "^.Ti (3.26) 
of the com6mec? jet-lepton distribution as follows: 



m.f± = . (3.27) 



Introducing the ratios 



^max 

rn ^ -J^ , (3.28) 

^max 

rf ^ > (3-29) 



we can now write the combined jet-lepton distributions for each lepton charge in terms of the 
unit-normalised variable ( p. 27] ) as 



\je+/s /=1J=1 



ifb, ifc) rl r'-^^p] {rlw?.i,+ ; x, y, z) 



2 2 ^ ^ 



7=1 J=l 



2 2 



dm^ _ i 2 



ji- J s 1=1 J=l 

^22 ^ 
+ 2^^^"^ ^'P^^'Pb.^c)r)j's^/]{r)m^j^^-,x,y,z). (3.31) 
/=i J=i 

Note that whenever the two endpoints and m^"^ are different, one of the two ratios 

r„ and is guaranteed to exceed 1, so that there will be a range of masses for which the 
corresponding argument (r„mj£ or rf-rhj^) in the J'g.^jj functions would exceed 1 as well. This 

is why it was necessary to extend the range of definition of our ^g.jj and ^g./j functions in 
Appendix]^ to be < m < oo, although it seems trivial, since the functions vanish identically 
for m > 1. 
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As can be readily seen from eqs. ( |3.22| ) and ( |3.23| ), both of these observable distributions 
are unit normalised 



dN 



dm 



dN 



dm?- 



(3.32) 
(3.33) 



just like the observable dilepton distribution ( p.?] ) (see eq. ( 2.44| )). 

This concludes the derivation of our first main result. It is worth recapitulating what 
we managed to achieve so far. We obtained exact analytical expressions for the three ex- 
perimentally observable invariant mass distributions: dilepton (2/7), jet plus positive lepton 
( 3.30 ) and jet plus negative lepton ( 3.31 ). All three of our formulas are unit normalised and 
can be readily rescaled for the actual observed number of events (which is the same for each 
of the three distributions). Our formulas are written in terms of a set of known functions 



Tgij which are explicitly defined in Appendix The coefficients Kjj appearing in our 
formulas are defined in eqs. ( gl6| - ^l9|) , (l3^p^ and (^37| - ^) , and depend on only three 
model-dependent parameters (pa, and ipc- Those parameters are defined in eqs. ( |1.9| ) and 
dSll^ll , and are a priori unknown, so that they must be measured from the data. 

The basic idea of our spin measurement method (whose main steps will be presented in 
detail in the next section) will be to fit our formulas to the shapes of the measured invariant 
mass distributions. Since there are 6 possible spin configurations, this fit will have to be 
repeated 6 times - once for each value of S. Since we have only three parametric degrees 
of freedom ipa, fb and (pc, with which we are trying to fit three whole distributions, one 
would expect that the fit will be successful only for the correct spin configuration S and for 
the remaining 5 spin cases the fit will fail. Indeed we find that this expectation is generally 
correct, and in Sec. || we shall give explicit examples of how this procedure might work in 
practice. However, we also find that there are two pairs of "twin" spin scenarios, discussed 
in Sec. 4.1, which are often completely indistinguishable, even as a matter of principle. 



3.2 Invariant mass formulas in the {J^sa^-^sp^-^s-^^-^s-s^ basis 

While the fitting exercise just described can in principle be performed with our results writ- 
ten in terms of the basis functions from Appendix we find that for the actual 
practical application of our method, it is much more convenient to rewrite our results in 
a different functional basis. We therefore introduce an alternative set of basis functions 
{•^s'i' -^5^3' •^s'-l' -^s^l} which are nothing but linear combinations of those appearing in our 
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old set: 



ip) 



t(p) _ 7r(P) I 7r(P) _ TTU 



ip) 



S;f3 



ip) 



J ( 



5:5 



■Trip) I -r-(p) _ -r-(p) _ 

"I" ■^S;12 ■^S;21 ■^S;22 ( ' 



-r-(P) _ •r-(P) _ -r-(P) i -r(P) \ 

-^5;12 -^5;21 -^5;22 f ' 



^ I S;ll "I" 5;12 "T 5;21 "T -^5;22 f ' 



r(p) 



r(p) 



r(p) 



(3.34) 
(3.35) 
(3.36) 
(3.37) 



for any p € {££,j£n,jif}- Using the normalisation conditions ( 2.10| - 2?T^ ), it is easy to see 



that the newly defined functions .^5?^, •^g'^ and J'g'l^ are zero- normalised 



J^^S^fjiiTi^ X, y, z) drh^ = , 
J^^ghrn^; x, y, z) dw? = , 



(p) 

while the function J-Q.i is unit-normalised 



> 

T'f\{rh?] X, y, z) drh^ = 1 . 



(3.38) 
(3.39) 
(3.40) 

(3.41) 



The explicit form of the new basis functions {-^s^L -^S' -^S' -^s!!} 

can be easily obtained 

by substituting the results from Appendix |S| into the definitions ( p. 34 - 3.371 ) . The result is 
given in Appendix ^. 

The advantage of the new set of basis functions becomes immediately apparent when we 
rewrite our results for the different invariant mass distributions: 



dN 
dmi, 



dN 



= (^a! 2;, y, z) + a{ipb, <fa)^s-l ^' ^) ' 



dw? 



i^n / S 



(3.42) 
(3.43) 



dlSl 
drh?. 



/ / s 



- J :Fg^s\w?j^^. ■ X, y, z) + a{(fb, fa)J='s-^J^ {m]f,^ ; x, y, z) 



±P{(pc,^b)J^sJ i'rnjef;x,y,z) ±'y{ipa,(pc)^s-J 2;, y, z) ^ , (3.44) 



where a, (3 and 7 are constant coefficients related to the chirality parameters (1.9) as follows 

Oi{^b,fa) = cos 299fc cos 2v9a , (3.45) 
Pi<fc,V^b) = cos 2(^c cos 2(^fc , (3.46) 
7iv^a,^c) = cos 299a COS 2(^c • (3.47) 
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Each one of the a, (5 and 7 parameters can take values in the interval [—1, 1]. However, a, (3 
and 7 are not completely unrelated. Given their definitions (3.45-3.47), it is easy to see that 
they must satisfy certain relations among themselves, and those are listed in Appendix |y. 

Using the normalisation conditions ( 3.38^3^411 ), one can easily show that all distributions 
(|3. 42^3.441 ) are properly normahsed as in eqs. ( ^.44 ^.22| , ^.23|) . Eqs. (|3.42| - |3.47l ) represent 
our main theoretical result. In the remainder of this section we shall discuss and interpret 
those equations. In the subsequent sections we shall illustrate how Eqs. ( 3.42| - 3^47| ) can be 
used for measurements of the spins, couplings and mixing angles. 

There are several desirable features of the {^'s-^a^-^^s-'p^-^'s-^^-^^s-^s^ basis used to write 

down eqs. ( ^.42 - 3.44 ). First, consider the •^5^] terms which appear without any parametric 

(p) 

coefficients. In most cases for S and p, the function J'g.^ simply gives the invariant mass 
distribution as predicted by pure phase space, i.e. where any spin correlations are ignored. 
This is true whenever there are only scalars and/or fermions among the intermediate particles 
appearing between the SM fermion pair whose invariant mass is being calculated. However, 
if a heavy vector boson appears among the intermediate heavy particles, the J'g.^ function 
always deviates from pure phase space. In fact this deviation cannot be compensated by a 
judicious choice of the a, (3 and 7 parameters. Therefore, one of our general conclusions will 
be that a heavy vector boson always leads to deviations from pure phase space and conversely, 
whenever a pure phase space distribution is observed, a heavy vector boson can be ruled out. 



Another nice feature of eqs. ( 3.42 - 3.44 ) is that the three parametric degrees of freedom are 
now explicit in terms of the coefficients a, (3 and 7. Even more importantly, it is immediately 
apparent which particular combination of the model-dependent parameters 99^5 and (pc 
(i.e. which combination of couplings and mixing angles) can be measured from any given 
distribution. For example, the observable dilepton invariant mass distribution L~g~ given in 
eq. ( 3.42| ) only depends on a, but does not depend on /3 and 7. Since the dilepton distribution 
is experimentally observable, this would allow a direct measurement of the a parameter from 
the dilepton data alone, by fitting to the shape predicted by (^.42)) . Note that a{(pi,,ipa) 
depends only on the chirality parameters ipi, and ipa entering the corresponding vertices for the 
near (^„) and far (£j) leptons. The fact that a (and as a consequence, the dilepton invariant 
mass shape ( ^.42 )) does not depend on the chirality parameter (pc associated with the quark 
vertex, should be intuitively obvious - the two leptons are not affected by the preceding 
events higher up in the cascade decay chain (see Fig. ||). The resulting measurement of a 
can be immediately interpreted in terms of the underlying chirality parameters ipa and ipb, 
as illustrated in Fig. |3|, leading to one constraint among ipa and 92;,. Clearly, the a{ipi),Lpa) 
measurement alone is not sufficient to pin down the precise values of ipa and (pi,. However, once 
it is supplemented with the additional measurements of (3{(pc-,^Pb) and j{ipa,<fc) as explained 
below, in principle all three parameters ipa, ft and ipc will be completely determined. 

Similarly, we can see that the jet-near lepton invariant mass distribution ( |3.43| ) only 
depends on the parameter /?, and does not contain the parameters a or 7. Again notice from 
Fig. 2 that f3{(pc,^b) in turn depends only on the chirality parameters (pc and ipb associated 
with the corresponding vertices for the quark (g) and the near lepton (£«)• This is also 
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Figure 3: The parameter a{(pb, ipa) defined in (3.45) as a function of cosipa and cosifb 



intuitively clear - the jet and near lepton should not be affected by what happens later down 
the decay chain. A measurement of (3 therefore can be immediately interpreted in terms of 
the underlying chirality parameters ipc and ifb, and the relationship is exactly the same as 
the one exhibited in Fig. ^ between a{(pb, fa) and its arguments. 

However, as we already explained in Sec. the {jin} invariant mass distribution ( |3.43| ) 
is not separately observable, and instead has to be combined with the {jif} distribution given 



in ( 3.44 ) to form the experimentally observable {ji } and distributions. We see from 

eq. ( 3.44| ) that the {j^f} distribution depends on all three parameters a, (3 and 7, which is 
again easy to understand intuitively - the intermediate lepton does affect its neighbors on 
both sides (g and if). Given the expressions ( ^l3D and (|3l4| ), we can immediately combine 
them using the same procedure as in eqs. (|3^ and (|3^ ): 



dN \ 



^ V n irnrri%± ; X, y, z) + rj J'^l^'' {r)w?.^± ; x, y, z) 



± -^sV (^/^i£±; y' ^) T .F^g^Km^V; ^> y. ^) • (3-48) 



Notice that the same j3 and 7 terms in ( 3.4g| ) appear with opposite signs in the {i^^} and 
the distribution. This suggests that instead of the two individual distributions ( |3.4g| ) 
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we should be considering their sum 

S^-{m^,;x,y,z,a)^\^\ +1^1 (3-49) 



clN \ ( dN 



2 -r-(i^n)/ 2-2 \ I 2 T-ii^f) t 2-2 \ , 2 -r-U^f) r 2-2 \ 

'^n^S;5 {r^mj^;x,y,z) + rjj'g./ [r frrijf.x.y, z) + ar f J'g.J {rfmjf,x,y, z) 



and their difference 

I), - K,; y> A 7) - j ^ - j ^ (3.50) 

= jrjj^^g.!l\rjm^jf;x,y,z) +Prjj^^g.^p^\rjm^jf;x,y,z) - Pr^J^^^.^p\rlm^jf,x,y,z) . 
The normahsation conditions for the newly defined quantities S^" and are 

Sg^ {m'^f ]x,y, z,a) dmj^ = 2 , (3.51) 

Dg^ {rh'jf,x,y, z, P,'~f) dm'ji = . (3.52) 
Eq. (3.49) reveals one of our most important results - that the sum of the two jet-lepton 







distributions depends on a single model-dependent parameter, and more importantly, this 
is the same parameter (a) which also determines the dilepton invariant mass distribution. 
Therefore, once a is measured from the relatively clean dilepton data, the experimentally 
observable Sg^ distribution is completely specified! This is a very important result, and as 
we shall see later in our examples, the dilepton [L^ ) and Sg^ distributions by themselves 
can often discriminate among the various spin alternatives. 

Of course, the Dg~ distribution is also observable, and it can be used as an additional 
cross-check of the results obtained with the two a-dependent distributions. The importance of 
the distribution is that it can provide a measurement of the other two model-dependent 
parameters f5 and 7. Note, however, that the 7 parameter can be measured only if 5 = 4, 5, 6, 
since for the remaining three cases we have 

4-7"^ = -^S'^ = for 5 = 1,2,3, 

and Dg^ becomes 7-independent. Similarly, the parameter /3 can only be determined for 
S = 1, 4, 5, 6 since for the remaining two cases 5 = 2,3 

4%"^=4%'^ = for 5 = 2, 3, 

and Dg~ becomes /3-independent as well. 

Now we are in a position to contrast our approach to previous spin discrimination studies 
based on the lepton charge asymmetry [19]. The latter is simply the ratio 

.+-.-2 n , Dt~{m]^;x,y,z,(3,j) 

{mji,;x,y,z,a,(3,-y) = , ' — . (3.53) 

6 J {m^-^;x,y,z,a) 
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We can immediately see that, in general, A'^~ is a much more model-dependent quantity than 
either Sg~ or . Indeed, as we just discussed, depends on a single model-dependent 
parameter (a), Dg~ depends on two other model-dependent parameters (/3 and 7), while, as 
evidenced by eq. ( |3.53| ), Ag~ depends on all three of these (a, /3 and 7). Second, the lepton 
charge asymmetry is not normalised to any particular constant numerical value, unlike the 
S'^^ and Dg~ distributions (see eqs. ( p.51 , 33^ )). But most importantly, Ag~ is a single 
distribution, derived from Sg~ and Dg~ , therefore it is bound to contain less information 
than the two separate distributions Sg~ and . Our explicit examples in Sec. |5| will show 
that, as might be expected, the useful information contained in is approximately the 
same as the information contained in . Therefore, by considering in addition the S^^ 
distribution, as we are suggesting here, one is recovering the information which was lost when 
forming the ratio (3.53). This information gain is most striking for the case of a pp collider 
like the Tevatron, as discussed in detail below in Sec. 4.2. 



4. The method 



The starting point in our analysis is the set of analytical formulas ( |3.42 , 3.49 , 3.50 ) derived 
in the previous section for the three experimentally observable invariant mass distributions: 
dilepton , and sum {Sg~) and difference {D'^~) of the and the distributions: 



{rhu;x,y,z,a) 
rO'^n)/ 2^2 .. 



dN 



dN 



i^ii'^ x,y,z) +a J^^*:^^ {mjf, x, y, z) , 



drh?:., j 
/ s 



+ 



dN 
drh?.,,. 



(4.1) 
(4.2) 



ri^s'i"' {i-lmjf^]x,y, z) + rjj'j^.^/\rjm]f,x,y,z) + arj J'^g^'J' {rjmjf,x,y, z) 



2 Tr0'^/)r^2 - 2^., 



Ds {rh'j(,;x,y,z,^,j) = 



dN \ 



dN \ 



(4.3) 



= 'y'^j^s-n\f}^'jeiX,y,z) + Prjj^g.^p^\rjmjf,x,y,z) -l5rlj^^^^p\rlT^^;x,y,z) . 

The functions T^^^, J^^*^, T^^^ and J^^^J are given in Appendix while the constant model- 
dependent parameters a, (3 and 7 were defined in eqs. ( 3.45 - 3.47 ): 

OL{^b,^a) = cos2ipbCos2ipa , (4.4) 
Pi'fc, ^b) = COS 2(^c COS 2ipb = if - f) COS 2(pc cos 2ipb , (4.5) 

lifa, 'fc) = COS 2ipa COS 2(fc = {f - f) COS 2ipa COS 2ipc , (4.6) 

where in the last two equations we have used the relation ( 3.15| ). The angles (fa, ^b and 
were defined in eq. ( [L.9| ) and parameterise the relative chirality of the corresponding inter- 
action vertex in Fig. while the particle-antiparticle fractions / and / were introduced in 
Sec. 1.2 and satisfy eq. ( |1.3D . Given the data for the three distributions (4.1-4^), one then 



- 28 - 



tries to fit for the unknown model-dependent coefficients a, /? and 7, considering eacli of the 
six different spin possibilities S one at a time. The result will be 6 different sets of "best fit" 
values for these coefficients, {as, Ps^ls}, S = {1,...,6}, and an accompanying measure for 
the goodness of fit in each case. The fits can be done simultaneously for all three parame- 
ters, or alternatively, one can first determine a from the relatively cleaner L~^~ sample, and 



subsequently use this fitted value of as in eqs. (4.2,4.3). The goodness of fit for each S will 
indicate whether this particular spin configuration is consistent with the data or not, and, 
given the expected experimental statistical and systematic errors, one can also readily assign 
confidence level probabilities to those statements. As we have been emphasizing throughout, 
this procedure is completely model- independent, and in fact produces an independent mea- 
surement of the model-dependent parameters a, (3 and 7, which can then be translated into a 
measurement of the underlying theoretical model parameters (pa, '■Pb, and /. For example, 
when all three parameters a, (3 and 7 are measured and found to be non-zero, one can invert 
eqs. (4.4-4^) and solve for (pa-, (pb and Lpc up to a two-fold ambiguity: 



COS2y?a = ib-^ V^a/57, 
1 



cos2(/Jb = ±- A/a/37 , 
7 

cos 2(pc = ±- J - \/a(ff , 

f - f Oi 



(4.7) 
(4.8) 
(4.9) 



where in all three equations one should take either the "+" or the "— " sign on the right-hand 
side. The origin of this two-fold ambiguity is easy to understand. Observe that the defining 



equations (4.4-4.6) for a, (3 and 7 are invariant under the simultaneous transformations 



vr 



92a , 



TT 



TT 



(4.10) 



whose effect is precisely to flip the signs in the right-hand sides of eqs. (4.7-19|). Given the 



defining relation ( |1.9D , the transformations (4.10) are equivalent to the chirality exchange 



IolI ^ IorI , I&lI ^ \bR\ , \cl\ 



(4.11) 



The physical meaning of eq. ( 4.11 ) is clear - we can only measure the chirality of the three 
different vertices in Fig. || only relative to each other. When choosing the plus signs in 
eqs. (|4.7| - ^ ), we get a solution for the couplings with one particular set of chiralities, while 



choosing the minus sign in eqs. (4^-0) yields a solution where the couplings have just the 
opposite chiralities. Since there is nothing to provide a reference point for the chiralities, it 
is impossible to remove this L ^ R ambiguity without making some model assumptions, or 
without considering additional independent measurements. Using the solutions ( [4.7[]4.9| ) and 
the definitions (|1.9[) we can write down the general solution for the couplings in terms of the 
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measured parameters a, (3 and 7, as 



\hR\ 
\cr\ 
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1 
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72 



1 



1 =F - \/a^ 
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1 1 , N 2 



/-/« 



(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 



where the appearance of the it sign is due to the two-fold ambiguity just discussed. Here 
the two solutions are obtained by choosing the upper or lower sign in each equation, corre- 
spondingly. It is worth making a few comments regarding eqs. ( [4.12 - 1.17 ), which represent 
our second main result. 



Note that while in general a, /3 and 7 can have either sign, eqs. ( 3.45 - 3.47 ) imply that 
the product a/37 is always non-negative. Furthermore, from eqs. ( p5H3l7| ) it also follows 
that \al3\ < I7I, |/?7| < \a\ and \'ya\ < Therefore all square roots in eqs. ( 4.12| - 4.17 ) 
are well behaved and never yield any imaginary solutions. It is interesting to note the 
dependence on the particle-antiparticle fraction / discussed in Sec. We see that for any 
given measurement of a, (3 and 7, the effective couplings |ol|, \aR\, I^lI and associated 
with the particle A and particle B vertices of Fig. |^ can be uniquely determined, up to the 
two- fold L ^ R ambiguity ( 4.111 ). In other words, the particle-antiparticle ambiguity T2 
discussed in the Introduction only affects the determination of the \cl\ and \cr\ couplings, 
as seen from eqs. ([4.16 - 07|) . The values of the couplings \cl\ and \cr\ are not uniquely 
determined, and instead are parameterised as a function of /. Although we do not know the 



exact value of /, consistency of eqs. ( [4.16H4.l7[ ) restricts the allowed values of / to be in the 
range 



< / < - 




1- \ — 



or 




(4.18) 



The fact that the allowed range for / splits into two separate intervals could already be 
seen in Fig. ||: notice that there are two disjoint branches in the (cosipcf) plane which are 
consistent with a given fixed value of (^ci i-e. with a given set of measured a, /? and 7. At a 
pp collider like the LHC, in general we expect / > so we should select the higher / range 
in eq. ( [4.18 ), while the lower / range in eq. ( 4.18[ ) would be relevant for a hypothetical pp 
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collider ("anti-LHC"): 



2 \ V a 



LHCM: + </<l, (4.19) 



anti - LHC (pp) : < / < ^ ^1 - ^J^j ■ (4.20) 

While eq. ( [4.19| ) is not a real measurement of the value of / at the LHC, it nevertheless 
contains very important information. For example, if the measured values of a, [3 and 7 
happen to be such that \f3'^\ « |a|, then / becomes very severely constrained, and the 
restriction ( ^.19| ) by itself is sufficient to yield a measurement of the value of /: / ~ 1. 

In the following Section |5| we shall give numerous examples of how our method might 
work in practice. But before we conclude this section we shall anticipate some general results 
which can be gleaned from our analytical formulas (^-4.3). In particular, in Sec. 4.1 



we 



shall show that the two pairs of spin configurations FSFS and FSFV, as well as FVFS and 
FVFV, very often give identical results for the invariant mass distributions, and cannot be 



differentiated without additional model assumptions. Then in Sec. [4.2| we shall show that our 
method is also applicable at the Tevatron, where in contrast the lepton charge asymmetry 
A'^~ is identically zero for all spin configurations S and thus contains no useful information. 

4.1 The twin spin scenarios FSFS/FSFV and FVFS/FVFV 

Consulting the definitions of the functions in Appendix one can see that 

= , (4.21) 

= = > (4-22) 
•^g = .Fg = , (4.23) 

•^g = -^g (4.24) 
for any p G {i£,j£n,jif}. Therefore the relation 

l-2z , , 

is sufficient to guarantee that all invariant mass distributions (4J-0) are exactly the same 
in the case oi S = 2 (FSFS) and 5 = 3 (FSFV): 

Lt' (^rnjf,x,y,z,a3^-^-^^ = {mjf,x,y,z,a3) , (4.26) 

x,y,z, «3 Y^p^) = {^"ji^ x,y,z,a3) , (4.27) 
{rh]f, x,y,z, (32,12) = D^' {rh]f, x,y,z,/33, 73) . (4.28) 
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Note that this exact duphcation occurs irrespective of the values of the other two model- 
dependent parameters (3 and 7. In other words, relations (4.26- 0^ ) hold identically for 
any values of the five parameters a^, P3, 73, P2 and 72. As long as eq. ( [4.25 ) is true, the 
FSFS and FSFV models will yield identical invariant mass distributions for , and 
. This observation has very important implications for the eventual outcome of the 
spin measurement, if the data happens to come from one of those models, since the exact 
duplication ( 4.26 - 4^2^ ) then threatens to jeopardize our ability to discriminate among them. 
However, as we shall now see, whether discrimination is possible or not, depends on the 
actual values of a and z. Recall that the a parameter is defined in the range [—1, 1], while 
z is defined in (0,1), and therefore so is the ratio | }^2z l- Then, for any given value of 
03 G [—1,1], Q2 as given by ( 4.25 ) falls into its definition window, and an exact duplication 



takes place. However, the reverse is not true: not every value of 02 would lead to a valid 
solution for 03 according to eq. ( [4.25| ), since for large enough values of |a2|, the value of |a3| 
would exceed 1, which is not allowed. 

Our conclusion therefore is that the issue of confusing the two models FSFS and FSFV 
depends on whether the data comes from FSFV and we are trying to fit it with FSFS, or 
whether the data comes from FSFS and we are trying to fit it with FSFV. In the former 
case the two models will always be confused with each other, while in the latter case, the 
confusion arises only if 02 happens to satisfy 



|a2| < 



1 - 2z 



l + 2z 



(4.29) 



A close inspection of Appendix ^ also reveals a similar problem with the FVFS and 
FVFV spin configurations (5 = 4 and S = 5). In this case, we notice the following relations 



Trip) 



-rip) 1 



5;/3 



5;7 



2z 



l + 2z 



■Trip) 
-^4;/3 



1 - 2z 



^'^ l + 2z 



-pip) 



(4.30) 
(4.31) 
(4.32) 
(4.33) 



for any p £ {ii,jin,jif}. Therefore, the relations 

1 - 2z 
l + 2z ' 



q;4 = 05 



Pa = /35 



74 = 75 



l-2z 
l + 2z 



(4.34) 
(4.35) 
(4.36) 



would once again guarantee that all invariant mass distributions (4.1-4^) are exactly the 
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same in these two cases: 



5^ 



l-2z 
l + 2z 

1 - 2z 
l + 2z 



m^jf,;x,y,z,P5,-f5j 



1 - 2z 



+ 2z 



Lt {'mje;x,y,z,a5) , 



St {m'jf,x,y,z,a5) 



Dt {m^ji;x,y,z, (35,^5) 



(4.37) 
(4.38) 
(4.39) 



Following the same logic as before, we conclude that whenever the data comes from FVFV, 
the model will always be confused with FVFS. However, if the data comes from FVFS, the 
confusion arises only if and 74 happen to satisfy 



|«4| < 



1 - 2z 



ml 



< 



l + 2z 
l-2z 



l + 2z 



(4.40) 



(4.41) 



In addition to these two equations, the values of 04, P4 and 74 must also satisfy the domain 
constraints ( C.2| - C.5| ) from Appendix |y. 

4.2 Spin determination at the Tevatron 

At a pp collider such as the Tevatron, the symmetry of the initial state implies 



(4.42) 



On the surface, it may appear that this constraint eliminates only one out of the four model- 
dependent degrees of freedom (/, and (pc) that we originally started with. However, 
as can be deduced from eqs. (|3.13|J3.14 ) and also seen from Fig. |2|, the constraint ( [4.42| ) in 
fact completely fixes the (pc parameter 



vr 

and result both /? and 7 vanish identically: 

/3 = 7 = 



In that case from eq. (4.3) we have 







and a similar result holds for the lepton charge asymmetry ( 3.53| ) 







(4.43) 

(4.44) 
(4.45) 

(4.46) 



We see that at the Tevatron we do not learn anything from either or from the lepton 

charge asymmetry Ag~ . However, our results for L'^^ and S^^ still hold, and contain non- 
trivial spin information, so that the spin analysis following our method can still be performed. 
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In fact, our method can already be tested in the top quark semileptonic and dilepton samples 
at the Tevatron by looking at the invariant mass distribution of the 6-jet and the lepton [53]. 
Indeed, our decay chain from Fig. |^ can be applied to top quark decays, for example by 
identifying C = t, B = and A = I'l, and reinterpreting in as the 6-jet and if as the lepton 
coming from W decay. In that case, the rribi distribution should be described by our formula 
( |02D for L^'. Alternatively, one can identify the particles in Fig. |l| as I? = t, C = W~^, 
B = Vi, q = b and in = i- In this case, the rribe distribution will be described by our formula 



( 3.43 ) applied for S" = 4 or 5 = 5. In any case, one should observe the characteristic m term 

r-(«) 
6:<5 



in the invariant mass distribution (see the definition of J'i^^^ in Table ^ or the definition of 



J'l g and -^5^ in Table ^j, which would signal that the W is spin 1 and therefore the top 
quark and the neutrino are both spin 1/2. 



5. Determination of spins and couplings: examples 

In this section we shall give an explicit demonstration how to apply our method in practice 
at the LHC. We shall work out in detail 6 different examples, namely, we shall assume in turn 
that the observed data is coming from each one of the six spin configurations from Table |l[ 
Then we shall ask the question whether this data is consistent with one of the remaining 5 
alternatives. 

Since we do not yet have real data available, we will have to use simulated data. We shall 
therefore have to pick some values for the mass spectrum, couplings and particle-antiparticle 
fraction, namely we shall have to fix the values of x, y, z, (fa, <Pb, '■Pc, and /. In order to 
allow comparisons to previous studies in the literature, we shall use the parameters of the 
SPSla study point in supersymmetry. However, as advertised, we shall still perform the spin 
measurements in a model-independent way, i.e. as soon as we simulate our "data", we shall 
immediately "forget" how it was generated, and shall treat it as coming from a "black box" 
such as the actual collider experiment. 

For the SPSla mass spectrum we take the values used in Refs. [20,21] 

ruA = 96 GeV, = 143 GeV, mc = 177 GeV, mz) = 537 GeV , (5.1) 

which translate into 

X = 0.109, y = 0.653, z = 0.451 . (5.2) 
SPSla is characterised by the following approximate values for the coupling constants 

aL = 0, aR = 1, hi = 0, hR = 1, cl = 1, cr = 0, (5.3) 

and particle-antiparticle fractions / and / at the LHC 



/ = 0.7, 7=0.3 



(5.4) 



The spectrum ( ^.l] ) results in the following kinematic endpoints^ 



m^'^^ = mi5Vx(l-y)(l-z) = 77.31 GeV , (5.5) 



mT^ = 771^5^(1 = 298.77 GeV , (5.6) 



max 



m'-rr = m£,v'(l -x)(l - z) = 375.76 GeV , (5.7) 



m^f = mDVi^-x)il-yz) = 425.94 GeV . (5.8) 

Since we assume that the spectrum has been measured, the values of these endpoints are 
also known in advance of the spin measurement. We are therefore still allowed to write the 
measured invariant mass distributions ( ^4.1^4.31 ) in terms of the dimensionless invariant masses 

(HI)- 

Substituting the SPSla parameter choice (|5.3| ) and (|5.4|) into the definitions ( ^.45D -( ^747|) 
yields the following values of our model-dependent parameters a, (3 and 7 

a = l, /3 = -0.4, 7 = -0.4. (5.9) 

Note that a = 1 necessarily implies /3 = 7, in accordance with eqs. ( |3.45| )-( p^47| ). 

Eq. (^j^) defines the input values of the model-dependent parameters used in our study. 
We should reiterate that there is nothing special about the SPSla parameter choice, and we 
could have used any other study point instead. 

Using our method, we shall now perform 6 different exercises of spin determination. For 
each exercise, we shall take the input "data" to be given in turn by one of the six models from 
Table ||. We shall then try to fit the "data" to each of the remaining 5 spin configurations, 
using our general analytical expressions (4.1-4^ with floating, a priori unknown, parameters 



a, (3 and 7. Although the fit can be done simultaneously for all three parameters a, [3 and 7, 
we shall perform it sequentially, using the fact that the L'^~ and S~^~ distributions depend 
only on a and not on /? and 7. Therefore, we shall start with the cleaner L'g^ sample and 
first determine the value of a, which we shall then use to compare the thus predicted S^^ 
distribution to the "data". Quite often, it will be already at this stage that one could rule 
out all but the correct spin configuration. We shall encounter such examples below as well. 
Sometimes, however, there may still be several alternatives left, in which case we need to 
also consider the -D"' distribution, where we fit for the values of the coefficients (3 and 7. 
Details of our fitting procedure and examples of some fits are presented in Appendix |C|. Our 
results are summarised in Figs. ^ ^ and ^, which show our results for the ^5 , and 
distributions, correspondingly. In each of Figs. ^, |5| and |6| the solid (magenta) lines in each 
panel represent the input invariant mass distribution {Lg~, Sg~ or Dg~ , as appropriate) 
from our simulated "data", for each of the 6 spin configurations: a) SFSF; b) FSFS; c) FSFV; 
d) FVFS; e) FVFV; f) SFVF. The other (dotted or dashed) hues are our best fits to this 
data, for each of the remaining 5 spin configurations from Table |l|. The color code is the 
following. If the trial model fits the input data perfectly, we use a dashed (green) line. If the 



^The kinematic endpoint mjil^ is only needed for the extraction of the mass spectrum, while the actual 
distribution is not needed for our study. 
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fit fails to match the input data, we use (color-coded) dotted lines. The best fit values of a, 
(3 and 7 for each case are also shown, except for those cases (labelled by "NA") where they 
are left undetermined by the fit. Dotted lines of the same color imply that they are identical 
to each other, yet different from the input "data". 

5.1 SFSF example {S = 1) 

For the SPSla parameters ( |5.2[j5.4| ) (or alternatively, (^)), eqs. ( |4.1[ - [4.3D predict the following 
observable invariant mass distributions for the SFSF model: 



St 



1 



f 2.810 
1.228 
^ -2.880 logm^ 



D1 



-0.668 + 2.002 m' 



-0.035 



m'^jf, < 0.632 
0.632 < < 0.653 
0.653 < , 

< 0.632 
0.632 < < 0.653 



(5.10) 
(5.11) 



(5.12) 



6.633 - 6.633 + 5.481 log 0.653 < 



These distributions are shown with solid (magenta) lines in Figs. |^(a), ^(a) and ^(a), re- 
spectively. Following our procedure described above, we first try to fit the dilepton data in 
Fig. §(a). Due to the presence of an intermediate scalar particle B, the distribution for 
the SFSF chain (S=l), is completely flat. However, that does not necessarily mean that the 
spin of particle B is determined to be zero. In fact, as seen from Fig. ^(a), all other spin 
configurations except for S = 6 (SFVF) can also fit this flat distribution, simply by choosing 
a vanishing a parameter. Even the case of S = 6 (SFVF), whose "best fit" prediction is 
different from the input data, may still be difficult to discriminate in practice, once we factor 
in the finite statistics, detector resolution and combinatorial backgrounds. The bad news, 
therefore, is that we cannot immediately determine the spins from the L"' distribution alone, 
but the good news is that, as anticipated, we got a measurement of the a parameter, which 
represents some combination of heavy particle couplings and mixing angles. 

At this point it is worth comparing our Fig. 1(a) to Fig. 2a in Ref. [21], where a very 
similar exercise was performed^''. The two results are quite different, for example we find 
that 4 out of the 5 "wrong" models can perfectly fit the dilepton "data", while in Ref. [21] 
all 6 models give distinct dilepton shapes. Of course, neither of the two results is wrong, and 
the difference simply arises due to our different philosophy. In Ref. [21] the parameters a, 
(3 and 7 (in our notation) were all kept fixed to the SPSla values ( |5.9| ), while here we are 
allowing them to float, since they would not have been measured in advance independently. 
As a result, we tend to get much more similar distributions, indicating that once we factor 



^Fig. 2a of Ref. [21] is simply the collection of all six solid (magenta) lines in our Fig. kla)-(f), i.e. our 



input "data" for the six different spin configurations, using the same fixed SPSla values (5.2 
dependent parameters. 




) for the model 
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Figure 4: Dilepton invariant mass distributions (Lg ). The solid (magenta) line in each plot repre- 
sents the input dilepton distribution from our simulated "data", for each of the 6 spin configurations: 
a) SFSF; b) FSFS; c) FSFV; d) FVFS; e) FVFV; f) SFVF. The other (dotted or dashed) lines are our 
best Gts to this data, for each of the remaining 5 spin configurations from Table The color code is 
the following. If the trial model fits the input data perfectly, we use a dashed (green ) line. If the fit 
fails to match the input data, we use (color-coded) dotted lines. The best fit value of a for each case 
is also shown, except for cases where it is left undetermined (NA). 
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Figure 5: The same as in Fig. but for instead of 



in the experimental realism, the actual spin measurements might be even more challenging 
than previously anticipated. 

Having extracted all the relevant information out of the L"* distribution, we now move 
on to studying the distribution. As we already explained in Sec. 3.2, the advantage of 
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Figure 6: The same as in Fig. |^ but for D+ instead of 



considering as opposed to each one of the individual distributions and is 

that only depends on exactly the same parameter a as the dilepton distribution L"' 
(see eq. (|4.2|)). Since we have just measured a by fitting to the L"' data, at this stage there 
are no free parameters left in the S'^ distribution, and it is uniquely predicted for each of 
the 5 "wrong" spin scenarios. In Fig. ^a) we plot the resulting predictions for the six spin 
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models, using in each case the corresponding value of a, which had been measured in the 
previous step from the L"' distribution. We see that the distribution can now further 
differentiate between different spin cases, e.g. it can rule out (in principle) the 5 = 4 and 
S = 5 (FVFS and FVFV) models. Interestingly, now S = 6 (SFVF) gives a perfect match, 
but fortunately, it has already been eliminated from consideration by the analysis of the L"' 
data at the previous step. Unfortunately, the "wrong" spin scenarios S = 2 and S = 3 (FSFS 
and FSFV) once again give a perfect match to the data, so that even after considering both 
L"' and 5"' , we are still left with 3 distinct possibilities for the spins of the heavy partners. 
As we shall see later from the other 5 exercises, the SFSF input "data" is somewhat of an 
unlucky case, since we end up with several spin models which perfectly fit both the L"' and 
data. More often than not, and by themselves should be sufficient to narrow 
down the spin configuration alternatives to a single one (or at most two, due to the "twin" 



spin scenarios discussed in Sec. 4.1) 



We are therefore forced to consider our third piece of data, the distribution (|4.3| ). 
This distribution does not depend on the previously fitted parameter a, and instead needs to 
be fitted with the other two model-dependent parameters, /3 and 7. Even though D'^ itself 
does not explicitly depend on a, the fit is nevertheless impacted by the measured value of a, 
as the latter determines the allowed range of values for /3 and 7 (see Appendix |C| for details). 
The results from our fitting exercise to the SFSF "data" are shown in Fig. ^(a). We 

see that can now eliminate the remaining two "wrong" spin scenarios S = 2 and S = 3 
(FSFS and FSFV) and as a result of all three types of fits, we are able to determine uniquely 
the spin chain as being S = 1 (SFSF). In addition, we were also able to obtain a measurement 
of the parameter /3, which carries information about the couplings and mixing angles of the 
heavy partners D, C and B. Unfortunately, the parameters a and 7 are not experimentally 
accessible in this case (5 = 1), since their corresponding basis functions J'l^}, and are 
identically zero for any p G j^n, j(^f} (see Appendix 

Having investigated both the and distributions, we do not need to consider the 
lepton charge asymmetry A'^ , which is simply the ratio of and 5"' (see eq. ( p.53| )). 
Numerically the asymmetry A'^ and the difference show a very similar pattern of their 
distributions, and thus provide roughly the same amount of information. However, as we 
emphasized in Sees. |3.2| and there are cases where the asymmetry A'^ (as well as ) 
does not play any role at all. The cases of S=2 (FSFS) and S=3 (FSFV) discussed in the 
next subsection actually provide such an example. 

In the remainder of this section, we shall repeat the exercise that we just went through, 
each time taking our "data" from a different spin configuration, and trying to fit to it the 
remaining^^ 5 spin possibilities. 

5.2 FSFS and FSFV examples (5 = 2,3) 

With the SPSla parameters ( |5.9D , eqs. (|]^^^ predict the following observable invariant 



'Obviously the "correct" spin configuration will always give a good fit to its own "data" 
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mass distributions for the FSFS model 



Ll 



^2 



■ 2m« 
( 2.898 

1.316 



m]^ < 0.632 
0.632 < < 0.653 



-16.583 + 16.583 - 16.583 log m% 0.653 < m^^ 



, 



(5.13) 
(5.14) 

(5.15) 



which are shown by the solid magenta lines in Figs. |^(b), |5|(b), and ^(b), correspondingly. 
Similarly, for the FSFV model we get 



si 



1.052 - 0.104 
r 2.815 

1.233 



m^^ < 0.632 
0.632 < < 0.653 



0.860 + 0.860 m^^ - 3.590 log 0.653 < m^^ , 



it 



, 



(5.16) 
(5.17) 

(5.18) 



which are shown with solid magenta lines in Figs. §(c), |5|(c), and ^(c), correspondingly. 
The distributions ( 5.15 - 51151 ) and ( 5.16 - 5^T^ ) will be the input sets of data for our next two 
exercises. 

Perhaps the most striking feature in each of the data sets ( ^.13 -5.15) and ( ^.16 -5.18) is 
that the distribution, and consequently, the lepton charge asymmetry , are both 

identically zero. Therefore, they do not convey any information about the spins, since any 
spin configuration can fit those distributions with the proper choice of parameters as shown 
in Figs, ^(b) and ^(c). This being the case, we should concentrate on the L"' and 
distributions. 

First we shall discuss the case when the data comes from the FSFS {S = 2) model, 
eqs. (|5t|-|t|). a gain, we begin our analysis with , which in this case shows very good 
discrimination, and can already rule out all of the "wrong" spin combinations. As explained 
in Sec. 4T, FSFS {S = 2) can sometimes be faked by the FSFV (S=3) model, but this could 



only happen if the a parameter in the data satisfies eq. (4.29), i.e. 



\0i2\ < 



1 - 2z 



l + 2z 



0.05 



(5.19) 



Since for SPSla a = 1 (see eq. ( |5.9D ), this condition is not satisfied and the FSFV model 
cannot fake the FSFS "data". This is confirmed by our result in Fig. 0(b). 

Since the L"' distribution alone already singles out the correct spin configuration, we 
do not even need to consider the distribution. It is worth pointing out, however, that 
in this ideal case also can rule out all "wrong" spin models, although the differences 
are not so pronounced as for L"' , and in reality are likely to be washed out. In summary. 
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the FSFS "data" can be unambiguously interpreted in relation to the spin issue, and we can 

also get a measurement of the parameter a. On the other hand, the parameters (3 and 7 will 

(p) (p) 

remain undetermined, since their corresponding basis functions ^^'^ •^2-7 identically 
zero for any p £ {M, jin, jif} (see Appendix^). 

Now we shall discuss the case when the data comes from the FSFV {S = 3) model, 
eqs. This will provide our first example where our spin measurement ends up 



being inconclusive, yielding two different, equally plausible, possibilities for the spin chain. 
This result should have already been anticipated, based on our general discussion in Sec. 
There we showed that for any given FSFV data, the FSFS model {S = 2) can always provide 
a perfect fit, and furthermore, the value of 02 that would be measured for the "twin" FSFS 
model is 

02 = as ~ 0.05 , (5.20) 

where we used the SPSla values for 03 = 1 and z = 0.451. Our numerical study explicitly 
confirms this general expectation as shown in Figs. |^(c), ^(c) and |6|(c). In addition, we 
checked that the mju distributions for those two "twin" spin models are also identical. 



5.3 FVFS and FVFV examples {S = 4, 5) 



In this subsection we discuss the case of the other "twin" spin pair from Sec. 4.1, namely 
5 = 4 and 5* = 5 (FVFS and FVFV). Using the SPSla values for the model-dependent 
parameters, we obtain the following distributions for the FVFS case 



Si 



0.492 + 1.016 m|^, 
( 2.307 + 3.455 m^^ 

1.028 + 0.577 m^^ 



4.553 



-42.563 - 12.368 + 54.931 mj^ 
7.871 + 90.785 m|^) log m^^ 



< 0.632 
0.632 < rh^.g < 0.653 



0.653 < , 



(5.21) 



(5.22) 



-0.22 + 0.616 
-0.092 + 0.212 
-3.087 + 3.087 



m^jf, < 0.632 
0.632 < < 0.653 



0.874 + 2.678 ) log m^^ 0.653 < 



(5.23) 



which are shown with the solid (magenta) lines in Figs. H(d), |5|(d), and P(d), correspondingly. 
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For the FVFV case we get 



0.974 + 0.053 
2.496 + 2.908 m2 

1.217 + 0.030 m2 



4.553 mj^ 



27.809 - 43.679 mj^ + 15.870 

+ ('l4.382 - 4.710 rfi^^) logm^^ 



-0.139 + 0.415 

-0.011 + 0.011 rfi^^ 

1.109 - 1.109 m^^ 

+ (^1.004 - 0.139 m^^) logm^^ 



< 0.632 
0.632 < m^.f, < 0.653 



0.653 < , 



m^-^ < 0.632 
0.632 < < 0.653 



0.653 < m^^ , 



(5.24) 
(5.25) 



(5.26) 



and those are shown by the solid (magenta) hnes in Figs. |^(e), ^(e) and ^(e). 

The end result of the two exercises is very similar to what we obtained in the previous 
subsection for the other "twin" model pair (FSFS and FSFV). It could have also been an- 
ticipated from our general discussion in Sec. |4.l| . When going in the forward direction, i.e. 
starting with the FVFS "data" and fitting to it the other 5 models, we do not encounter 
any spin ambiguities. As already determined in the previous subsection, this is because the 
SPSla value of the a parameter (a = 1) does not satisfy the necessary condition ( [4.4C| ) for 
an FVFV model to fake the FVFS data. As a result, the two and distributions are 
already sufficient to pin down the spin case scenario, and the D'^ distribution can then be 
used as a cross-check and for a measurement of the (5 and 7 parameters. 

However, when going in the reverse direction, i.e. starting with the FVFV "data" and 
fitting the other 5 models including FVFS to it, we do encounter a spin ambiguity, just like 
in the S = ?> exercise above. Again, the reason for this was already explained in Sec. 4.1. In 
agreement with our analytical results, Figs, ^e), |5|(e) and ^(e) show that the FVFS model 
provides an identical match to the FVFV "data" for all three observable distributions L"' , 
and . The good news, however, is that while we are left with a two-fold ambiguity 
with respect to the spins, for each spin scenario the parameters a, (3 and 7 are precisely 
measured, so that we have independent measurements of three different combinations of the 



heavy partner couplings and mixing angles. In Sec. p.5| below we shall show how to interpret 
those measurements in terms of the more fundamental model parameters a^, a/j, 6l, hn, cl, 
CR and /. 



5.4 SFVF example (S = 6) 

Our final example is the SFVF spin chain, for which the SPSla model parameters (|57 
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predict the following observable distributions 



Li 



"^6 



0.405 TTi^^ , 



1.626 - 0.981 rfi^ 
{ 2.87 

1.288 



-0.344 - 4.493 + 4.837 



(5.27) 

m^^ < 0.632 
0.632 < m^-i, < 0.653 (5.28) 
5.870 log rfi^^ 0.653 < m|^, 



-0.322 



0.786 



-0.406 + 1.051 



m^jf, < 0.632 
0.632 < < 0.653 



5.870 

+ ( 3.384 



11.674m2^ + 5.804 



(5.29) 



3.595 



log m% 



J lug '"'ji 0.653 < iTT'jij 

shown with solid (magenta) lines in Figs. |^(f), |5|(f) and ^f). 

The case of S = 6 (SFVF) is very special, since in this case the dilepton invariant mass 
distribution ( ^.271) exhibits a characteristic term which is not present for any of the other 
5 spin configurations that we are considering. Note that the existence of an m'^ term in the 
dilepton SFVF data is generic, i.e. does not depend on the values of the model-dependent 
parameters such as a. This could be easily understood by realising that the dependence 

(££) 

originates from the "phase space" basis function J^^.g , which enters our general formula 
( |4.1| ) for the dilepton distribution without any model-dependent coefficients. More generally, 
an inspection of Table ^ reveals that the dilepton invariant mass distribution is in general 
given by some polynomial in terms of in^, whose power is equal to twice the spin of the 
intermediate particle B}^ Only in the SFVF case {S=6) do we have a spin 1 intermediate 
particle which brings about an term in . If the presence of this term can be observed 
in the dilepton data, it would unambiguously^^ signal the presence of a spin 1 mediator. Of 
course, the size of the coefficient of the term depends on the mass spectrum in the model, 
but it cannot be vanishingly small - this would require either z = 1 or y = 1, which would 
correspondingly close off the B ^ Ai oi the C Bi decay, and the whole decay chain will 
become unobservable. One of our general conclusions, therefore, is that the SFVF model^^, 
if it exists, should be discernible from the dilepton data alone. Our numerical results in 
Fig. Q(f) confirm this conclusion - we see that none of the other five models can reproduce 
the SFVF dilepton data, due to the presence of the term. Pictorially this can be seen 
from the fact that the predictions of the S = 1,2,3,4,5 models in Fig. |^ are always 
straight lines, while for the S = 6 model (SFVF) the prediction is never a straight line, due 
to the higher order m dependence. 

Before we move on to the next subsection, where we shall interpret our measurements of 
the a, P and 7 parameters, we briefly summarise the results from the preceding six exercises 



^■^A similar statement can be made about the niji^ invariant mass distributions from Table relating the 
power of the rh^e^ to the spin of the intermediate C particle. 

^^This observation is subject to our assumption that we do not consider heavy particles of spin 3/2 or 
higher. In general, an dependence would imply that the spin of the mediating particle is at least 1. 

^*0r more generally, a spin 1 or higher intermediate particle. 
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Result from 




Data originating from model 




Total No. 
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FVFV 
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L"^ only 
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1 


1 


4 


1 


2 


1 


4 


L+- © S+- © D+- 


1 


1 


2 


1 


2 


1 


2 



Table 3: Summary of the results from our spin discrimination analysis. Each entry represents the 
total number of models n which can perfectly fit the data sets listed in the first column, i.e. each 
entry n implies an n-fold model ambiguity of the corresponding data. The last column lists the total 
number of "wrong" spin configurations allowed by the corresponding data set, which was obtained by 
summing all the n's from the preceding 6 columns and subtracting 6 for the correct configurations. 



in Table |3|. The table shows the number of different spin configurations from Table |^ which 
can fit perfectly a given data set {L^ , , , or some combination thereof). Since this 
number depends on the spin configuration of the input "data" , we show 6 different columns, 
one for each different spin configuration of the "data" . The last column lists the total number 
of "wrong" spin configurations allowed by the corresponding data set in all 6 exercises. This 
number was obtained simply by summing all the entries from the preceding 6 columns and 
subtracting 6 to exclude the correct configurations among them. 

While one should be mindful that the number counts exhibited in Table ^ are only 
valid for the SPSla parameter choice, there are still some interesting conclusions which can 
be drawn from it. For example, we do not notice any particular pattern in the horizontal 
direction. In particular, the discriminating power of the different data sets, say L"' , and 

, varies greatly from model to model. There are cases where a single distribution works 
very well, for example L+- for FSFS and SFVF, S+- for FSFS and FVFS and D+- for SFSF, 
FVFS and SFVF. In all those cases the spin configuration is uniquely fixed by studying a 
single distribution! On the other hand, there are also cases where each one of these individual 
distributions performs rather poorly, for example L"' for SFSF, for SFSF and FSFV 
and for FSFS and FSFV. In the end, each one of the , and distributions, 
when considered in isolation, yields on the order of 10 fake spin configurations. What this 
simply means is that no single distribution can be universally "better" than the others. 

Things begin to get more interesting when we start combining information gained from 

2 or more different distributions. For example, when we combine any 2 out of our three 
observable distributions L"' , and , the total number of fake solutions drops down to 

3 or 4. Now again, which particular pair works better, is a model-dependent issue: ©S^ 
fails for the SFSF model, while L+- © D+- and 5+" © D+- both fail for the FSFV model. 
Finally, combining the information from all three distributions, L"' (BS^ (BD^ , we narrow 
down the remaining spin choices even further, but as we saw in Sees. 5.2 and 5.3, there are still 
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two cases of exact duplication, which are nothing but the "twin" spin scenarios of Sec. 4.1. 
Since this duplication is due to an exact mathematical identity, it will obviously still persist if 
we were to repeat our analysis including all the experimental realism (backgrounds, resolution, 
combinatorics, etc.). In fact, due to the expected imperfections in the real data, one may get 
even more duplicate examples, if anything. 
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Table 4: Available measurements of the model-dependent 
parameters a, (3 and 7 for each of the six spin configuarions. 



5.5 Measurements of couplings and mixing angles 

Recall that our general method 
from Sec. ^yields not only a determi- 
nation of a possible spin chain fitting 
the data, but also a measurement of 
the model-dependent a, (3 and 7 pa- 
rameters from eqs. (U^). Even m 
the spin duplication scenarios found 
in Sees. ^]2| and |5.3| , we still have a 
certain measurement of the a, (3 and 
7 coefficients for each of the two al- 
lowed spin chains. This is illustrated 
in Table Q where we summarize the 
available measurements of the a, (3 and 7 parameters in each individual spin case. Notice 
that only in the last three spin cases (FVFS, FVFV and SFVF) we are able to measure the 
complete set of all three parameters a, (3 and 7. In contrast, for the SFSF model chain we 
can only determine (3, while a and 7 remain unknown. On the other hand, for the FSFS and 
FSFV chains we can only determine a, while (3 and 7 remain arbitrary. In the remainder of 
this section we shall discuss the interpretation of those measurements in terms of the cou- 
plings and mixing angles of the heavy partners, i.e. we shall relate the measured values of a 
and/or (3 and/or 7 to the underlying model parameters /, ifa, and (pc- 

First, let us consider the case where we determine a spin chain to be one of the following 
three: FVFS, FVFV or SFVF. Then, as seen from Table ^ we will be able to measure the 
values of all three parameters a, (3 and 7. If we have correctly determined the spin chain, these 
values will be simply the starting SPSla inputs (|5.9|). Substituting those in eqs. ( [4.12 - 4.17 ), 
we obtain the two sets of solutions discussed at the beginning of Section ^ 



IolI = 0, \aR\ = 1, = 0, \hR\ = 1, \cl\ 



1 

2 + 2/ 



0.2 



\cr\ 



0.2 

2/-1 



,(5.30) 



and 



IolI = 1, \aR\ = 0, = 1, \hR\ = 0, \cl\ 



0.2 



2/-1 



\cr\ 



1 0.2 

2 + 2/^ 



-,(5.31) 



where the first (second) solution corresponds to choosing the upper (lower) sign in eqs. ( [4.12 - 
4.17] ). As expected, we obtain that each set is a one-parameter family of solutions, parame- 
terised by the value of the particle- antiparticle ratio /. The first solution set ( 5.30| ) reproduces 
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the SPSla parameter set for / = 0.7, but of course, we would have no way of knowing that 
/ = 0.7 is the correct value of /, since we would have to measure / independently by some 
other means. However, notice that even though we do not know the exact value of / at this 
point, the solutions ( 5.3C| - 5.31 ) unambiguously restrict the allowed range for / from ( [4.19| ) to 
be 

0.7 < / < 1 , (5.32) 

which is by itself already an important and useful experimental determination. 

Now let us discuss more specifically the case where the data is due to an FVFS or an 
FVFV spin chain (S* = 4 or S" = 5). As already explained in Sec. 4.1 and explicitly seen in 
our examples in Sec. fi.3[ here we may encounter a second solution for the spin chain, with its 
own measured a, (3 and 7 parameters. We remind the reader that when the data comes from 
an FVFV chain, there is always a duplicate spin solution due to an FVFS chain, while if the 
data comes from an FVFS chain, the duplicate FVFV solution exists only if the conditions 



( 4.40 , 4.41 ) are satisfied. While the duplicate spin chain prevents us from uniquely resolving 
the spin question, the interpretation of its a, (3 and 7 parameters can be done in a very 
similar fashion. Consider our duplication example from Sec. 5^ where an FVFS {S = 4) spin 
chain was able to "fake" the FVFV {S = 5) data. All three parameters a, f3 and 7 were still 
uniquely measured but the obtained values were not the starting SPSla values. Instead, our 
fitting procedure found 

a = 0.05, (3 = -0.4, 7 = -0.02 (5.33) 

as shown in Fig. ^(e) , H(e) and ^(e) . We see that the (3 parameter for the twin spin chain was 
found to be the same as the true (3 parameter of the data {(3 = —0.4), while both the a and 
7 parameters of the twin spin case are a factor of 20 smaller than the original inputs ( |5.9| ). 
This fact can be easily understood from our general results from Sec. 4.1 - the conditions 



(4.40, 4.41) which guarantee the existence of a duplicate solution, relate the values of a and 7 
for the two spin chains, and by the same factor of j^ff ~ where we have used the SPSla 
value of 2; = 0.451. Just as before, the measurements (|5.33|) translate into a measurement 



of the effective couplings and mixing angles as a function of /, up to a two- fold ambiguity. 
Substituting (|5.33|) in eqs. (4.12-4.17), we obtain the two solutions 



Iflil = 0.69, \aR\ = 0.72, |6l| = 0, \bR\ = 1, \cl\ 



1 0.2 



1 0.2 



2 2/-1' 
(5.34) 



or 



\aL\ = 0.72, \aR\ = 0.69, = 1, = 0, \cl\ 



0.2 



2/-r 



\cr\ 



1 

2 + 2/ 



0.2 



- 1 

(5.35) 

As expected, these solutions exhibit the same L <^ R symmetry ( [4.11 ) as the solutions (|5.30| ) 
and ( ^.31] ) for the "correct" spin configuration. Comparing eqs. ( |5.34| , 5.35| ) to eqs. ( 5.3C| , 
5.31j ), we see that we obtain the same result for the \cl\ and \cii\ couplings! In other 

words, although it may not be clear what is the correct spin chain - FVFS or FVFV, the 
chirality of the couplings at the quark and at the near lepton vertex will be known (up to the 
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inescapable two- fold ambiguity due to ( [4.11| )). This can be simply understood by noticing 
from eqs. ( 4.14 - 07| ) that the couplings \bji\, \cl\ and \cji\ only depend on a and 7 
through their ratio, which is the same for the correct and the fake spin solution, since a and 
7 are scaled by the same factor (see eqs. ( 4.34 , 06| ). Just as before, for the "wrong" spin 
chain we also obtain a constraint on the allowed range of the particle-antiparticle fraction / 
at the LHC: 

0.7 < / < 1 . (5.36) 



Notice that this is identical to the result ( p. 32 ) for the "correct" spin chain, so that the 
experimental determination of the range of the / parameter also does not suffer from the 
duplicate spin ambiguity. 

This concludes our discussion of the spin cases where we can measure all three parameters 
a, [3 and 7. For the remaining three spin chains, only partial information will be available 
(see Table |^. For example, in case of SFSF we can only measure the /3 parameter, which 
gives us one relation among ipi, and (pc 



cos2ipb cos2(^c = —0.4 , 
or alternatively, among and /: 

(2/ — 1) cos2(/7fe cos2(/3c = —0.4 . 



(5.37) 



(5.38) 



Unfortunately, we are unable to pin down further the precise values of (pc and /, and 
furthermore, ifa remains completely unknown. 

Similarly, in case of FSFS and FSFV, we can measure the a parameter, which gives us 
a relation between and 



a = cos 2(^6 cos 2ipa = 1 . 



(5.39) 



Normally, we would not be able to go any further, but the SPSla parameter set is "lucky" in 
the sense that it yields one of the two extreme values of a (see Fig. ^). In those circumstances, 
we can determine the actual values of (fa and 995, and subsequently, la^l, |aij|, |5l| and 
up to the usual L <^ R ambiguity: 



TT 



or 



\aL\=0, \aR\ = 1, =0, = 1 



la^l = 1, |a^j| = 0, = 1, =0 



(5.40) 
(5.41) 



Unfortunately, in either case, \cl\, |cr| and / will remain unconstrained. 

Finally, we briefly comment on the possibility of spin duplication between FSFS and 



FSFV discussed in Sec. 4.1 and Sec. 5.2. Here we will also obtain a measurement of the 



a parameter for the "wrong" spin chain. The two a parameters (for the "wrong" and for 
the "correct" spin configurations) are related according to ( 4.25| ) and the analysis for the 
couplings in the case of the "wrong" spin chain can be done in complete analogy. 
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6. Summary and conclusions 



We conclude by summarizing the main steps of our method for measuring spins, couplings 
and mixing angles of heavy partners in cascade decays with missing energy. We shall then 
contrast it to other proposals for spin measurements in the literature. 
The method involves the following basic steps. 

1. Data preparation. Identify a decay chain of interest which would yield three observable 
SM fermions. (In this paper we considered the example of a quark jet followed by 
two leptons, which is commonly encountered in models of supersymmetry and extra 
dimensions.) Then form the three observable invariant mass distributions for each pair 
of well-defined objects: {i~^£~}, and In order to remove the combinatorial 

ambiguities, perform an opposite- flavor subtraction on the leptons and a mixed-event 
subtraction on the jet. Apply final cuts to possibly suppress any SM and new physics 
backgrounds. As the end product from this step one obtains the three ditributions 
L+-, S+- and D+- defined in eqs. 



2. Mass measurements. This step is optional, since the mass measurements can in principle 
be performed simultaneously with the spin fits described below. However, in practice 
we expect that the invariant mass distributions would reveal their kinematic endpoints 
rather early on, so that the mass spectrum can be measured in advance of the spin 
determination. At the end of this step one would know the mass spectrum, i.e. the 
values oi X, y and z which enter the functions J^, as well as the kinematic endpoints 



^max -^j-^ig]^ yj^j^ normalise our invariant mass variables (see eqs. (|2.2| ) and ( p.27D ) 



3. Spin measurements. This step represents the actual spin measurement. One tries to 
fit^^ the data for the , and distributions obtained in Step 1 with the 

theoretical predictions (4.1-4^), for each value of S, i.e. for each set of allowed spin 



configurations for particles D, C, B and A (see Table If the fit is good, that 
particular spin chain is ruled in, while if the fit is bad, that particular spin chain will 
be ruled out. Our expectations for the generic outcome of this exercise are summarised 
in Table ^. When using the data from all three distributions , and , we 
expect that the fits will be able to rule out all but the correct spin configuration. The 



only exceptions are the spin duplication cases discussed in Sec. LI, when one may end 
up with at most two spin chain alternatives. 

4. Measurements of couplings and mixing angles. In this step one uses any available best- 
fit values for a, /3 and 7 obtained in the previous step, and determines the couplings 
{o-lI, l^i?!; I^lI; \t>B.\, and \cr\ from eqs. ( [4.12 - 4.17 ). There will be two different 



solutions due to the L ^ R symmetry, as discussed in Sec. |^ and illustrated with some 



^^In general, those are three-parameter fits for the floating, a priori unknown, coefficients q, /3 and 7. 
However, as discussed in Sec. ^ and iUustrated with our numerical examples in Sec. one could make use of 
the fact that the L"' and S'"' distributions only depend on the parameter a. Thus one could first extract a 
from and/or S^^, and then use this value to fit D^^ for /3 and 7, as shown in Appendix ^ 
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Data 




Can this data be fitted by model 




from 


SFSF 


FSFS 


FSFV 


FVFS 


FVFV 


SFVF 


OJj or 


yes 




Tl i~l 
LLKJ 






no 


FSFS 


no 


yes 


maybe 


no 


no 


no 


FSFV 


no 


yes 


yes 


no 


no 


no 


FVFS 


no 


no 


no 


yes 


maybe 


no 


FVFV 


no 


no 


no 


yes 


yes 


no 


SFVF 


no 


no 


no 


no 


no 


yes 



Table 5: Expected outcomes from our spin discrimination analysis, barring numerical accidents due 
to very special mass spectra. The two cases labelled "maybe" correspond to the potential confusion 
of an FSFS (FVFS) chain with an FSFV (FVFV) chain, which occurs only for a certain range of the 
model-dependent parameters - see eqs. ([l.29|) and (|4.40, f.4f ). 



examples in Sec. 5^. In addition, eq. ( 4.19| ) provides a restriction on the allowed range 



of values for the particle-antiparticle fraction / at the LHC. 

Having summarised the main steps of our method, we are ready to compare it to other 
approaches for spin measurements which already exist in the literature. In principle, no single 
method is universally applicable, therefore the availability of different and complementary 
techniques is an important virtue. Which method ends up being most successful in practice, 
will depend on the specific new physics scenario that we may encounter. With those caveats, 
we should point out some features of our method which are likely to make it relevant and 
successful, if a missing energy signal of new physics is seen at the LHC and/or the Tevatron. 

• Many of the existing techniques for spin determinations (see, for example, [30,31,45,46]) 
have been originally developed in the context of lepton colliders, where the total center 
of mass energy in each event is known. Consequently, at hadron colliders, those meth- 
ods are applicable only if the events can be fully reconstructed. In new physics scenarios 
with dark matter WIMPs, this appears to be rather challenging, since there are two in- 
visible WIMP particles escaping the detector. In some special circumstances, where two 
sufficiently long decay chains can be identified in the event, full reconstruction might be 
possible [10-12], but in any case, this appears to require very large data samples. In con- 
trast, our method relies on invariant mass distributions, which are frame- independent, 
and we do not need to have the event fully reconstructed. Furthermore, the event re- 
construction techniques currently being discussed rely on the pair-production of two 
heavy particles, both of which decay visibly to the lightest WIMP. Our method, on the 
other hand, does not require the presence of two separate decay chains in the event, 
and can be in principle also applied to the associated production of a WIMP with only 
one other heavy partner. 

• The invariant mass distributions L"* , and D'^ that we propose to analyse, are 
the basic starting point for any precision study of new physics parameters. In the past 
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they have been extensively discussed in relation to mass measurements, and we now 
simply propose to fully analyse them for the encoded spin information as well. 



One major advantage of our method in comparison to various event counting techniques 
[33,35,44,47] is that we do not need to know anything about a number of additional 
and a priori also unknown quantities such as the production cross-sections for the 
different parton-level initial states, the branching fractions, the experimental efficiencies, 
etc. Indeed, our method in essence only uses unit-normalised distributions, and is not 
affected by any of these additional variables. 

The previous three advantages are common to all studies which have relied exclusively 
on invariant mass distributions for spin determinations [19-22,32,36-43]. In comparison 
to those works, the main advantage of our approach is that it is completely general and 
model-independent, in particular we make no a priori assumptions about the type of 
couplings in each vertex of Fig. |l|, or about the particle-antiparticle fraction /. As a 
result, we were actually able to come up with measurements of certain combinations of 
those couplings and the / parameter (see Sees. ^ and 5.5). 



In conclusion, we reiterate that our goal in this paper was simply to present the basic 
idea of our method, and demonstrate that it can work as a matter of principle. Therefore in 
our analysis in Sec. ^ we did not include any realistic detector simulation, backgrounds (SM 
and combinatorial) etc. All of these factors will be investigated in a future publication [18]. 
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A. Appendix: The basis functions J^^^jj 

The basis functions J^g.fj {rh?]x, y, z) are listed in Table ^ and the basis functions T^^^^ j{m? ; x , y, z) 
are given in Table ^. Below we explicitly show the remaining basis functions J^g.^/j {m^; x, y, z) : 



SFSF {S = 1) 



J'^^.'(\m^;x,y,z) = r^^.[i\rh'-x,y, z) 



FSFS {S = 2) 



J^illi\m^ ■,x,y,z) = ^2-21 {rn^;x,y,z) 



^2^12 i^"^ ;x,y,z) = J^2;2^^ {m'^;x,y,z) 



(1-2/)^ 



' (1 — y) + log y if < y 
1 — + logm^ if y < rfi^ < 1 (A.l) 

if > 1 

'(l-y) + ylogy if < y 

1 — + 2/ log if y < < 1 (A. 2) 

if > 1 

' (1 — y) + log y if < y 

1 — + logm^ if y < rfi^ < 1 (A. 3) 

if > 1 

(1 — y) + y log y if tv? < y 

1 — + y log if y < < 1 (A. 4) 
if > 1 



FSFV {S = 3) 



•^3;ii ("^ ;x,y,z) =T^.2i ;x,y,z) 

' (l-y)(l-2z) + (l-2yz)logy if < y 
(1 - m^){l - 2z) + (1 - 2yz)logm2 if y < < i (A.5) 
if m2 > 1 



(l-y)2(l + 2z) 



•^3;12 ^ (^-^ ^' ^) = ^3-22 ' l"^"! y^ ^) 

' (1 — y)(l — 2z) + (y — 2z) log y liiv? <y 

(1 - m2)(l - 2z) + (y - 2z) log if y < < 1 (A.6) 

if > 1 



(l-y)2(l + 2z) 
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FVFS {S = 4) 



(l + 2x)(2+y)(l-y)2 



' {l~y)[4x-y- 4m2(2 - 3.x)] 
+ [(-1 + 4x)y + 4m2{l - {2 + y){l - x)}]logy if < y 

(1 - m'^)[Ax{2y + 1) - 5y - Am^{l - x)] 

+ [(-1 + 4a;)y + 4m2{l - (2 + y)(l - x)}]\ogw? if y < m? < 1 
if > 1 



(A.7) 



T'ilf{n?\x,y,z) 



6 



(l + 2x)(2 + y)(l-y)2 



' (1 - y)[2 + 3y - 2x(5 + y) + 4m2(2 - 3x)] 
+[y(4 + y) - 4a;(l + 2y) - A.Tr?{\ - (2 + y)(l - a;)}] log y lirr? <y 

(1 - m2)[2 + 9y - 2a;(5 + 6y) + 2^2(1 - x)\ (A.8) 
+[y(4 + y) - 4a;(l + 2y) - 4vr?{\ - (2 + y)(l - a;)}] logm^ if y < < 1 

if > 1 



T'i^-^{vr?\x,y,z) 



(l + 2x)(2+y)(l-y)2 



6 



(l + 2x)(2+y)(l-y)2 



(l-y)[-y-4m2(2-x)] 

— [y + 4?ft-^{l + y(l — x)}] logy if < y 

(1 - ir?)\-hy — 4m^(l - x)\ 

— [y + 4m^{l + y(l — x)}] logm^ \i y <rr? <\ 

if > 1 



' (1 - y)[2 + 3y + 2x(l - y) + 4m2(2 - x)] 
+[y(4 + y) + 4m^{l + y(l - x)}] log y if < y 

(1 - rr?)\l + 9y + 2x(l - 2y) + 2m2(l - x)] 

+[y(4 + y) + 4m^{l + y(l — x)}] log rr? if y < < 1 

if > 1 



(A.9) 



(A. 10) 



- 53 - 



FVFV {S = 5) 



J^^^.^(\m';x,y,z) - + 2a;)(2 + y)(l - + 2^) 



( (1 - y)[4:X - y + 2z{2 + 3y - 2x{5 + y)} 
-4m2(2 - 3x)(l - 2z)] - [y - 2yz{A + y) + 4x{2z - y{l - 4z)} 
+4m2{l +y-x{2 + y)}{l - 2z)] log y 



X < 



(1 - rn^)[Ax{l - 5z + 2y{l - 3z)} by + 2z{2 + 9y) 
-4m^{l - x){l- z)] - [y- 2yz(4 + y) + Ax{2z - y{l - Az)} 
+Am^{l + y - x{2 + y)}{l -22;)] log 







if m? < y 

(A.ll) 

if ?/ < TV? < 1 

if > 1 



6 



if rh < y 



X < 



(l + 2x)(2 + y)(l-2/)2(l + 2z) 

' (1 - ?/)[2 + 3y - 2x{5 + y) + 2{4x - y)z 
+Am?{2 - 3x)(l - 2z)] - [4x{l + 2y{l - z)} -y{4. + y- 2z) 
-Aw?{l + y-x{2 + y)}{l - 2z)] log y 

(1 - rh?)[2 - 2x{5 - Az + 2y{3 - iz)} + y{9 - lOz) 
+2m2(l - x){l - 4z)] - [4x{l + 2y{l - z)} -y{A + y- 2z) 
-4m^{l + y- x{2 + y)}(l - 2z)] log if y < < 1 

if > 1 

6 



(A.12) 



(l + 2x)(2 + y)(l-y)2(l + 2z) 

' (1 - y)[-y + 2{2 + 2x(l - y) + 3y}z - 4m^{2 - x){l - 2z)] 
-[y{l - 2(4 + y)2;} + Am?{l + y- xy){l - 2z)] logy \im? <y 

X < {l-w?)[4.{l+x)z-y{b-2{9-Ax)z}-Aw?{l-x){l- z)] (A.13) 
-[y{\ - 2{4. + y)z} + A!m^{l + y - xy){l - 2z)]\ogm^ if y < m? < 1 

if rh^ > 1 

6 



X < 



(l + 2x)(2 + y)(l-y)2(l + 2z) 

' (1 - y)[2 + 2x(l -y)+ y{3 - 2z) + \rr?{2 - x)(l - 2z)\ 
+ [y(4 -\-y — 2z)+ 4'm?{l + y — xy){l — 2z)] logy if < y 

(1 - m2)[2 + 2x{l - 2y) + y(9 - lOz) + 2m^{l - x){l - 4z)] (A.14) 
+ [y(4 + y - 2z) + Aw?{l + y - xy){l — 2z)] logw? ily <m? <\ 

if > 1 
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SFVF (5 = 6) 



(l+2y)(l-y)2(2 + z) 



6 

(l + 2y)(l-y)2(2 + z) 



(l+2y)(l-y)2(2 + z) 



6 

(l+2y)(l-y)2(2 + z) 



' (1 - y)[2 - 3z - 2y(l + z) + 4m2(l - 2z)] 

— [z(l + 4y) — 4777.2(1 — z — yz)] logy if < y 

(1 - m2)[2 -3z- 8yz + 2771^(1 - z)] 

— [z{l + 4y) — 4m2(l — z — yz)] logm^ if y < rn^ < 1 

if > 1 



(A.15) 



' (1 - y)[2 - 3z + 2y(5 - z) + Am'^{3 - 2z)] 
-[z{l + 4y) - 4y(2 + y) - 4rn2(l + 2y - z - yz)] logy ilm? <y 

(1 - m2)[2 - 3^ + 4y(5 - 2z) + 2rfi2(l - z)] (A.16) 
— [z(l + 4y) — 4y(2 + y) — 4rn2(l + 2y — z — yz)] \ogw? if y < < 1 

if > 1 



' (l-y)[^-4?7i2(l-2z)] 
+ [yz — 4m?{l — z — yz)] logy if m? < y 

(1 - m'^)[z{l + 4y) - 4m2(l - z)] 

+ [yz — 4?77^(1 — z — yz)] logm^ if y < w-^ < 1 

if > 1 



' (1 - y)[z - 4y - 4m2(3 - 2z)] 

— [y(4 — z) + 4r7T.2(l + 2y — z — yz)] logy ii < y 

(1 - m'^)[z - 4y(3 - - 4777,2(1 - z)] 

— [y(4 — 2;) + 4771^(1 + 2y — 2; — yz)] log777.2 if y < rn^ < 1 

if 771^ > 1 



(A.17) 



(A. 18) 
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s 


Spins 


^sfi (^'; ^. y. ^) = ^s;i2 (^'; ^. y. ^) 


•^5;2i^ (m2; a;, y, 2;) = J^g^^ (m^; y, z) 


1 


SFSF 


2m2 


2(1 - m^) 


2 


FSFS 


1 


1 


3 


FSFV 


1 


1 


4 


FVFS 


(1+2XK2+?/) iy + 4(1 - y + - 4(1 - .t)(1 - y)mi} 


(1+2^(2+^)^43^ + y + 4(1 - 2z - y + xy)m^ - 4(1 - a;)(l - y)m'^} 


5 


FVFV 


(1+2 j(2+y) {y + 4(1 - y + xy)m2 - 4(1 - x){l - y)m^} 


(i+2x){2+v) {4^^ + y + 4(1 - 2.T - y + xy)m2 - 4(1 - a;)(l - y)m'^} 


6 


SFVF 


^{2y+(l-2y)m2} 


- (1 - 2y)m2} 



Table 6: Basis functions for the j£n invariant mass distribution. 



s 


Spins 




f "2 \ tW ( -2 \ 


1 


SFSF 


1 


1 


2 


FSFS 


2(1 - rn^) 




3 


FSFV 


^{1 - (1 - 2z)m'} 


^{2^ + (l-22;)m2} 


4 


FVFS 


^{y + (2 - y)m2} 


^^{2 - (2 - y)mn 


5 


FVFV 


(2+y)ti+2z){y + + (2 - y)(l - 2z)m^} 


(2+,)a+2.){2 + 2y2 (2 2.)m2} 


6 


SFVF 


(i+2^)(2+.){4y + ^ + 4(1 - 2y - z + yz)m^ - 4(1 - y)(l - z)m^} 


(1+2^(2+^)1^ + 4(1 - z + yz)m2 - 4(1 - y)(l - z)m^} 



Table 7: Basis functions for the dilepton invariant mass distribution. 



B. Appendix: The basis functions T^gl^, ^s^p, ^s'-\ -^s^l 



The basis functions J'g^^, ^s^p^ -^5^7 -^^s^s listed in Table ^. The basis functions 
J^s-o^^ ' •^5'^"'*' -^s-s"^ given in Table ^. Below we explicitly show the remaining 

basis functions T^i^J'\ J^^J.^p''\ t''^^^'^ and J'^^^/^: 
SFSF (S = 1) 



J^[^.g^\'ifP;x,y,z) 



FSFS {S = 2) 



T^\'j\m^-x,y,z) = ^ 



yf 



' 2(1 + + logy if < y 

2(1 — m^) + (1 + ?/) log \i y <rh? <1 

if > 1 

' (1 — y) log y if < y 
(1 — y) log if y < "T-^ < 1 
if > 1 

' 2(1 - y) + (1 + y) log y if < y 

2(1 — m^) + (1 + y) log if y < rfi^ < 1 

if > 1 



r\^'f\m^;x,y,z) = T%'\rn'-x,y,z)=^ 



-3-0'^/)/ - 2 ^ 



(i-y)^ 



' (1 — y) log y if < y 
(1 — y) log if y < w-^ < 1 
if > 1 



FSFV (5 = 3) 

^^f;j^^(m2;x,y,z) 



-1 l-2z 
(l-y)2l + 2z 



(B.l) 
(B.2) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



' 2(1 -y) + (l + y) logy if < y 
2(1 - m2) + (1 +y)logm2 if y < < 1 (B.7) 
if > 1 



^g^^(m2;x,y,z) = J^^'f\m^;x,y,z)=0 

' (1 — y) log y if rfi^ < y 



(m ;x,y,z) 



(1 — y) log if y < "rh^ < 1 
if > 1 



(B.8) 



(B.9) 
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FVFS {S = 4) 



(l+2x)(2 + y)(l-y)2 

—6x 

(l+2x)(2 + y)(l-y)2 



6x 

(l + 2x)(2 + y)(l-y)2 

-7-0'^/)/ "2 \ 

J^l/'im ■,x,y,z) 
3 

(l + 2x)(2 + y)(l-2/)2 



' (1 - y)[-2(l + 2y) + 2x(3 + y) - 16m2(l - x)] 
-[y(5 + y)-2x(l + 3y) + 8m2(l-x)(l + 2/)]logy if < y 

(1 - m2)[-2(l + 7?/) + 6x(l + 2y) ~ 6m'^{l - x)] (B.IO) 
-[y(5 + y) - 2x(l + 3y) + 8m2(l - x){l + ?/)] logm^ if y < < 1 

if > 1 



f 2(1- y) + (l + y) logy if < y 
2(1 — m^) + (1 + y) log 777,2 if y < < 1 
if > 1 



f 2(l-y)(l+y)(l-x) 
+[-2x(l + y)+y(3 + y)]logy if < y 

2(1 - m2)(l - x){(l + 2y) - m^} 

+[-2x(l + y) + y(3 + y)] log if y < < 1 

if > 1 



(B.ll) 



' 4(1 - y)(l + rn^) + [(1 + 3y) + 4m^] log y iim^ <y 
4(1 - m2)(l + y) + [(1 + 3y) + 4m2] logm^ if y < < 1 (B.12) 
if m2 > 1 



(B.13) 



- 58 - 



FVFV (5 = 5) 



1 - 2z 



(l + 2x)(2 + y)(l-y)2l + 2z 



' (1 - y)[-2{l + 2y) + 2x(3 + y) - lQm?{l - x)\ 
-[y{b + y)-2x{l + 2,y) 

+8rn2(l — + J/)] logy if < y 

(1 - m?)[-2{l + 7y) + 6x(l + 2y) - 6m'^{l - x)] (B.14) 
-[y(5 + y)-2x(l + 3y) 

+8m^{l — x){l + y)] log rh? ii y < < 1 

if > 1 



J^l\^p^\w?\x,y,z) 



—Qx 

(l+2x)(2 + y)(l-y)2 



6a; 



' 2(1 -y) + (l + y) logy if < y 

2(1 - rh^) + (1 + y) log rfi^ if y < < 1 
^0 if > 1 



(l+2x)(2 + y)(l-y)2 l + 2z 



'4(l-y)(l+m2) 
+ [(1 + 3y) + 4m2] logy iiw?<y 

4(l-m2)(l+y) 

+ [(1 + 3y) + Aiv?] \ogrh? \iy <m? <l 
if > 1 



\^ ;x,y,z) 



(l+2x)(2 + y)(l-y)2 



f 2(l-y)(l + y)(l-x) 
+[-2x(l + y)+y(3 + y)]logy if < y 

2(1 - m2)(l - x){(l + 2y) - m^} 

+[-2x(l + y) + y(3 + y)] log if y < < 1 

if > 1 



(B.15) 



(B.16) 



(B.17) 
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SFVF (5 = 6) 





(l+2y)(l-y)2(2 + z) 



(l+2y)(l-y)2(2 + z) 



(l + 2y)(l-y)2(2 + ^) 

-7-0'^/)/ '2 \ 

J^Q./'ijn ■,x,y,z) 
3 

(l + 2y)(l-y)2(2 + z) 



2(1 -y) + (l + y) logy if < y 
2(1 - m^) + (1 + y) log if y < < 1 
if > 1 



(B.18) 



' (1 - y)[2(l + 3y) - 2(2 + y)z + 16m2(l - z)] 
+ [2y(3 + y) - (1 + 5y)z + 8m2(l + y)(l - z)] logy if < y 

(1 - m2)[2(l + 8y) - 4(1 + 3y)z + 6m'^{l - z)] (B.19) 
+ [2y(3 + y) - (1 + 5y)z + 8m2(l + y)(l - z)] logm? if y < rh^ < 1 

if > 1 



' 4(1 - y)(y + m^) + [y(3 + y) + 4ym^] logy if < y 
8y(l - m^) + [y(3 + y) + 4ym2] log ni^ if y < < 1 (B.20) 
if rh^ > 1 



f 2(l-y)(l + y)(l-2) 
+[-(1 - y)(l + 2y) + (1 + 3y)(l - z)] logy if < y 

2(1 - m^){l - z){{l + 2y) - m^} 

+ [-(1 - y)(l + 2y) + (1 + 3y)(l - z)] logm^ if y < < 1 
if > 1 



(B.21) 
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s 


S|")ius 








T^.. {nr-.x.ij. z) 


1 


SFSF 


1 











2 


FSFS 


1 


1 - 2m^ 








3 


FSFV 


1 










4 


FVFS 


1 


-1^^(1-2-') 








5 


FVFV 


1 


(2-y)(l-2z) o^2x 
{2+y)(l+2z)y^ ) 








6 


SFVF 


(i+2y)(2+z){^y + ^ + 4(1 - 2/)(l - z){m^ ~ m^)} 


il+2y)(2+z)^^ > 









Table 8: Basis functions for the dilepton invariant mass distribution. 



s 


Spins 






J-g"^(m2;x,2/,z) 


.Fg'^"^(m2;a;,y,z) 


1 


SFSF 


1 





-(1 - 2m2) 





2 


FSFS 


1 











3 


FSFV 


1 


(J 





(J 


4 


FVFS 


(i+2x)(2+i,){2a^ + y + 4(1 - x){l - y){m^ - m*} 





TiT2lfey(l 2rn^) 





5 


FVFV 


(l+2c^){2+y){'^^ + y + ^i^ 





(i+2x)(2+?;)(l 2m) 





6 


SFVF 


1 





-i^(l-2rn^) 






Table 9: Basis functions for the jin invariant mass distribution. 



C. Appendix: Fitting procedure for the parameters a, (5 and 7 

In the absence of any error bars, we use a rather naive matching criterion, namely 

X^(a,/3,7) = y (/o("^^ao,/3o,7o) - /("^^a,/3,7)^ dm^ , (C.l) 

where /o(m^, ao) /^O) 7o) represents the experimental data that needs to be fitted and /(m^, a, /3, 7) 
is the theoretical prediction for it. We then minimize the x^{a^ (3, 7) function for a, (5 and/or 
7, as appropriate, ao, /?o and 70 are fixed constant values of the a, j3 and 7 parameters 
as predicted for the corresponding study point. A more sophisticated analysis including the 
expected statistical uncertainties is postponed for a future publication [18]. 

As we discussed in Sec. ^ fitting to the or to the 5"' distribution is a simple one- 
parameter fit for a, while fitting to the data is a two-parameter fit for j3 and 7. Fig. |^ 
shows sample results from our fits for (3 and 7 performed in the course of the exercises 
described in Sec. ^. In each plot in Fig. |^ the "data" /o(m^, ao, /3o> 7o) comes from the first 
spin chain (shown in red) at the top of each plot, which is then fitted with the distribution 
/(m^, a, (3, 7) predicted by the second spin chain (shown in blue). The contour lines represent 
constant values of x^(q^i/3i7)) where a has already been fixed by fitting to . The blue 
dot corresponds to the absolute minimum of x^, ignoring any restrictions on a, (3 and 7. 
However, the parameters a, (3 and 7 are not completely independent from each other. For 
any given a, the physically allowed region in the (/3,7) parameter space is described by an 



envelope which satisfies 

a/3 < 7, /37 < a, 7a < /? , if a > 0, ;9 > and 7 > , (C.2) 

a/3 > 7, /37 < «, 7a > /? , if a > 0, /? < and 7 < , (C.3) 

a/3 > 7, /37 > a, 7a < /? , if a < 0, /3 > and 7 < , (C.4) 

a/3 < 7, /37 > Q, 7a > /? , if a < 0, /3 < and 7 > . (C.5) 



In Fig. ^ we denote this allowed region in white (sometimes it may reduce to a single line). 
The green triangle corresponds to the minimum of the function within this restricted 
parameter set. The green triangle solution for j3 and 7 was then used for our plots in Fig. ^ 
For the two cases with FVFV (S=5) "data", the global minimum happens to lie within the 
(white) allowed region and so the blue dot and the green triangle coincide. 

For the extreme values of |a|, the (white) allowed region collapses to one or two lines: 

/3 = 0or7 = 0, ifa = 0, (C.6) 
7 = ±/3, ifa = ±l. (C.7) 
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